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In a model of an Abelian vector boson with a Maxwell kinetic term and non- negative mass-squared 
it is demonstrated that, under fairly general conditions during inflation, a scale- invariant spectrum 
of perturbations for the components of a vector field, massive or not, whose kinetic function (and 
mass) is modulated by the inflaton field is an attractor solution. If the field is massless, or if 
it remains light until the end of inflation, this attractor solution also generates anisotropic stress, 
which can render inflation weakly anisotropic. The above two characteristics of the attractor solution 
can source (independently or combined together) signiflcant statistical anisotropy in the curvature 
perturbation, which may well be observable in the near future. 
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I. INTRODUCTION 



Cosmic inflation is arguably the most compelling way to overcome (or at least ameliorate) the so-called horizon and 
flatness fine-tuning problems for the initial conditions of the Hot Big Bang cosmology (H- However, this success of 
■ inflation is guaranteed if inflation lasts long enough regardless of the specifics which cause the accelerated expansion 
C^' in the first place. As a result one cannot use the horizon and fiatness problems as a discriminator between different 
^ ^ models of infiation. To this end, much attention has been paid to another byproduct of an inflationary period; namely 
1^ ' the generation of the curvature perturbation C in the Universe, which is the source of structure formation and the 
. observed CMB primordial temperature perturbation. The characteristics of ^ open an observational window on the 
dynamics of inflation, which is no more that ten e- folds and corresponds to the so-called cosmological scales 0. 
\ Using this, inflation has become a testable theory with each model realisation offering specific predictions that can 
J> ' be observationally falsified. 

t — \ Contrast with observations has strengthened the case for inflation. Firstly, in the last decade, the rival theory (that 
of cosmic strings) for the generation of C, was observationally rejected Q. Secondly, the basic predictions of inflation 
' appear to be supported by observational evidence. Indeed, inflation generically seems to produce a predominantly 
scale-invariant spectrum of predominantly Gaussian curvature perturbations which are also predominantly statistically 
homogeneous and isotropic Q. Simple, single- field infiationary models are able to produce such a C,. However, in 
recent years the precision of cosmological observations has increased to the point that allows us to explore deviations 
from the above "vanilla" predictions of inflation, which has put increasing tension on inflationary models and resulted 
in a shift in the thinking of model-builders, away from the simple single-fleld toy-models and more towards more 
complex model realisations of inflation which make better use of the rich content of realistic theories beyond the 
standard model. A prominent example is the so-called curvaton hypothesis which removes the responsibility for 
the generation of C, from the inflaton fleld and assigns it to another field (the curvaton) which is unrelated to the 
' physics of infiation (it does not affect the inflationary dynamics) and can be associated with physics at much lower 
energy scales, e.g. TeV physics explored by collider experiments such as the LHC. 

Observations now indicate that there is a significant deviation from scale-invariance in C,. Indeed, the spectral index 
of the curvature perturbation is ris — 1 = — 0.037 ± 0.014 according to the latest WMAP data [H. This observation 
reveals the existence of infiationary dynamics, which agrees with expectations. It also means that models of inflation 
that generate a blue spectrum of perturbations are most probably rejected. In conjunction with the upper bounds 
on tensor perturbations, the above result has also falsified some red-tilted models too, such as all monomial models 
of chaotic inflation higher than quadratic. Similarly, non-Gaussian signatures appear to offer a second test for 
inflationary models. Indeed, the latest WMAP results produce a tentative observation of non-Gaussianity in the 
squeezed configuration, which gives for the non-linearity parameter /nl = 32 ± 21 If this result is confirmed it 
will falsify all single-field models of infiation, which predict |/nl| ^ 1- 
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Along these lines deviations from the remaining "vanilla" predictions of inflation should be investigated, namely 
statistical homogeneity and isotropy. Observations already offer some tentative evidence that ^ may not be exactly 
homogeneous and isotropic. For example, a 10% difference in power between hemispheres has been reported which 
would amount to a deviation from homogeneity Q. Furthermore, tantalising evidence of a preferred direction (the 
so-called Axis of Evil 0) on the microwave sky, amounting to an alignment of the low multipoles of the CMB, has been 
much discussed Q • This alignment persists beyond foreground removal @ and is statistically extremely unlikely ■ 
If such a preferred direction does exist then this could amount to statistical anisotropy in the curvature perturbation, 
which can become a powerful discriminator between inflationary models and mechanisms for the generation of C. 

Observations currently allow up to 30% statistical anisotropy in the power-spectrum of C [llj.^ Similarly, statistical 
anisotropy in the bispectrum can be predominantly anisotropic even if statistical anisotropy in the spectrum is 
negligible. If this is the case, then non-Gaussianity will feature an angular dependence on the microwave sky [l3l [l^. 
The forthcoming observations of the Planck satellite may well provide conclusive evidence of non-zero statistical 
anisotropy in the spectrum and/or bispectrum of the CMB temperature perturbations. In particular, the bound on 
statistical anisotropy in the spectrum will be reduced to 2% if such anisotropy is not observed [l^. Similarly, the 
observational bounds on /nl are expected to tighten by at least an order of magnitude, so statistical anisotropy in 
the bispectrum may well be found if non-Gaussianity is observed. Thus, it is imperative to investigate whether and 
how inflation can generate a statistically anisotropic signal. 

Traditionally, the dynamics of inflation and the generation of C. have been studied using fundamental scalar fields 
only. Such a setup, however, cannot produce a preferred direction on the microwave sky. The simplest way to do so 
is to involve at least one vector field, which naturally selects a preferred direction. Using a vector field to contribute 
to the curvature perturbation was pioneered by Ref. [l^ , where it was shown that a massive Abelian vector field can 
successfully play the role of the curvaton. Since then, substantial interest was ignited for the involvement of vector 
fields in inflation and the generation of the curvature perturbation. The first comprehensive study of how vector fields 
can contribute to C and generate statistical anisotropy can be found in Ref. [l7| . 

Involving a vector field can give rise to statistically anisotropic ( in two ways. One possibility is that the contribution 
of the vector field to the density during inflation may result in non-negligible anisotropic stress. This can produce 
statistically anisotropic perturbations of the inflaton scalar field which can eventually give rise to statistical anisotropy 
in ( |18l - l20| . In this case, the perturbations of the vector field itself are unimportant (and can be ignored) and the 
only requirement is that its contribution to the energy density remains non-negligible throughout infiation, or at least 
during the period when the cosmological scales exit the horizon. The second possibility is orthogonal to the previous 
one, in that it assumes a negligible anisotropic stress during infiation and it focuses instead on the contribution to 
^ of the vector field perturbations, which can give rise to statistical anisotropy if the particle production process for 
the vector field is itself anisotropic [l3, [HI- This is the generic expectation for a vector field that undergoes particle 
production and obtains a superhorizon spectrum of perturbations for some or all of its components [T7| . 

To undergo particle production and obtain superhorizon spectra of perturbations for its components, the vector 
field needs to have its conformal invariance appropriately broken during inflation. Mechanisms for such a breaking of 
conformality were originally developed in the effort to generate a coherent primordial magnetic field during inflation 
(e.g. see Ref. j22l. [23| and references therein). One proposal suggested couplin g no n-minimally an Abelian vector field 
to gravity through a term of the form ^RA^, where R is the scalar curvature [2J|. The model was shown to generate 
scale-invariant spectra for all vector field components (introducing a small non-zero mass gives rise to a longitudinal 
component) and particle production was anisotropic at the level of 100% [l3, 125]. Thus, if such a vector field generates 
a subdominant contribution to C^M '^^^ generate observable statistical anisotropy. This model, however, may suffer 
from instabilities such as ghosts [2^ (see however Ref. (27|). 

To overcome the instability issue another model was put forward, which employs a vector field with a varying 
kinetic function (2l] |. This possibility has been massively explored for the generation of a primordial magnetic field 
(28| . In this model, particle production of the transverse components of the vector field can produce a scale-invariant 
superhorizon spectrum of perturbations if the kinetic function of the vector field scales as / cx a~^^'^ during inflation 
(or at least when the cosmological scales exit the horizon), where a is the scale factor of the Universe [2l|, |23. If the 
vector field has a non-zero mass m it makes sense to talk about the longitudinal component as well. To obtain a scale 
invariant spectrum for this one, the additional requirement is m oc a |2l| . We also need the mass of the physical (in 
contrast to comoving) vector field AI = m/^/J to be small (i.e. M H, with H being the Hubble parameter) when 
the cosmological scales exit the horizon. Under the above conditions, if the vector field remains light until the end of 



^ In Ref. Illl statistical anisotropy of magnitude (29 it 3)% is actually found in the spectrum of the WMAP data at 9-ct. However, the 
preferred direction is suspiciously close to the ecliptic plane so the authors acknowledge that this observation is probably influenced by 
some unknown systematic (see also Ref. Il2ll ). Hence, this finding can only be taken as an upper bound on the statistical anisotropy of 
the primordial signal. 
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inflation then the power spectra of the transverse and longitudinal components respectively are \2l\ 

H\ . „ f HY f H 



so that Va, ^ and particle production is strongly anisotropic. In contrast, if the vector field becomes heavy by the 
end of inflation (i.e. M > H) then the vector field begins coherent oscillations around zero and the power spectra of 
the transverse and longitudinal components obtain the average value [2l| 



which means that particle production becomes isotropic. This possibility arises only if AI is growing during inflation, 
which corresponds to / oc (the other possibility, namely f (x. a? , results in M = constant). One can translate the 
above spectra of vector field perturbations into a contribution to C, through a variety of mechanisms. In Ref. [2l| the 
vector curvaton mechanism of Ref. [3 used. In Ref. [lOl the end of inflation mechanism was used instead for the 
same vector field model with m = 0. 

One expects the modulation of / and m to be due to a degree of freedom which varies during infiation. The most 
natural candidate for this is the inflaton field itself. However, in this case, there is a coupling between the vector field 
and the inflaton, which can backreact onto the inflaton's dynamics. Such backreaction was first considered in Ref. (l9| 
where it was shown that, for a specific functional form of the kinetic function, a model of quadratic chaotic inflation 
renders the scaling / oc a^^ an attractor solution. The role of the backreaction to vector field particle production 
in quadratic chaotic inflation with a massless, Abelian gauge fleld was also investigated in Ref. [31| with the aim to 
explore the generation of a primordial magnetic fleld. The same attractor behaviour was found. Furthermore, the 
backreaction on inflation from a massless, Abelian vector field whose kinetic function is modulated by the infiaton 
was studied in Ref. (s^ . where an attractor behaviour which results in weakly anisotropic inflation was also observed. 

In this paper, we perform a stability analysis to obtain the criteria under which the conditions / oc and to oc a, 
which render the vector field perturbation spectra scale invariant (and which at first glance may seem somewhat 
artificial), become attractor solutions for a generic functional form of the dependence of the kinetic function f{(j}) and 
the scalar potential V{4)) on the infiaton field 0. We find that the conditions imposed on the model parameters are 
such that the above scalings for / and m become attractor solutions for 50% of the parameter space (in the sense 
that for a given model parameter the condition is one upper or lower bound only) and can be attained for natural 
values. The criterion boils down to the requirement that the backreaction onto the inflaton field dynamics should not 
be decreasing during inflation, because were it so, the vector field could never affect the inflationary dynamics and no 
attractor behaviour would arise. Thus, we show that the conditions which guarantee the generation of scale invariant 
spectra of vector fleld perturbations (with amplitudes as shown above) can be naturally attained if the kinetic function 
and mass of the vector field are modulated by the infiaton during infiation. 

Our findings also show that the attractor solution can generate a non-zero anisotropic stress when the vector field 
is still light. This anisotropic stress can be kept small (negligible) by choosing the model parameters appropriately. 
In this case, the assumption of an approximately isotropic infiationary expansion isiustified and we can rely on the 
vector field perturbations to produce statistical anisotropy as discussed in Ref. (TtI . [21| . However, we can choose 
the model parameters such that the anisotropic stress of the attractor solution is non-negligible in which case the 
inflationary expansion itself would be anisotropic. In this case, the vector fleld perturbations can play a negligible 
role (and can be ignored) and we can produce the statistical anisotropy from the anisotropic expansion itself as in 
Refs. [l^ [23|. In our paper we do not make a choice between these two options, or indeed the third possibility in 
which both the above two sources of statistical anisotropy can be important. Instead we quantify the anisotropic 
stress at the attractor solutions as well as demonstrate that the scaling of / and to results to scale invariance. 

An interesting byproduct of the vector field backreaction on the inflationary dynamics is that it provides extra 
resistance in the roll of the inflaton down its potential. This can be used to maintain slow-roll even if the curvature 
of the scalar potential along the inflaton direction is substantial. Thereby we can overcome the infamous ?7-problem 
of inflation which arises by Kahler corrections to the scalar potential in supergravity theories. 

The structure of our paper is as follows. In Sec. [H] we outline our model and obtain the equations of motion for 
all the fields involved. In Sec. IIIII we define a set of dimensionless quantities which we use to perform a dynamical 
analysis of our model, without specifying the functional dependence of /(</>), to(0) and V{(t)). In Sec. IIVI we study 
the case of a massless vector field, in which only the transverse components are physical. We obtain the standard 
slow- roll critical point, which is an attractor for half of the parameter space where the vector field backreaction does 
not grow during inflation. We also flnd another critical point which is our vector scaling solution, which exists for 
the remaining half of the parameter space, when the vector field backreaction to the inflaton does affect the inflaton's 
roll, while also rendering the standard slow-roll critical point unstable. In Sec.|V]we consider a non-zero mass for the 
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vector field and investigate whether our vector scahng solution is spoiled or not. Firstly, we study the case of a light 
vector field (M ^ H) where we find that our vector scaling solution trifurcates into three critical points with weak 
flows between them. The inflationary expansion at these points can be anisotropic. Afterwards, we study the case 
of a heavy vector field (M > H) which is undergoing harmonic oscillations. Again we find a vector scaling solution, 
which now corresponds to zero anisotropy. In Sec. IVII we display a summary of our results. In Sec. IVIl] we apply our 
findings to four different concrete examples and we numerically verify our analytic results. The first example assumes 
an exponential dependence of /(</)), m(0) and V{(f>) on the inflaton field, which could be motivated by string theory 
considerations. The second example considers chaotic infiation and reproduces the findings of Ref. jlj in the massless 
field case. The third and fourth examples are supergravity inspired and assume that the infiaton is a flat direction 
of supersymmetry, lifted only by Kahler corrections to the scalar potential. The vector field is then a Higgsed gauge 
field with / being the gauge kinetic function associated with the gauge coupling as / ~ 1/5^- In these examples we 
clearly demonstrate how the backreaction of the vector field can overcome the yy-problem of infiation in supergravity. 
Finally, in Sec. IVIIll we discuss our results and present our conclusions. 

While this paper was being written up Ref. [33 appeared, which performs a similar stability analysis for the massless 
version of our model assuming an exponential functional dependence for /((/>) and V^(0). As such it corresponds to 
part of the first of our examples. However, the authors of Ref. (33| have explored deeper the inflationary dynamics and 
were mostly interested in the possibility of anisotropic inflationary expansion, completely ignoring the perturbations 
of the vector field and what the vector scaling solution implies for them. Our results agree with their findings where 
there is overlap. 

In our paper we consider natural units, where c = h ~ = 1 and Newton's gravitational constant is SttG = J^ip^, 
with mp being the reduced Planck mass. 



II. EQUATIONS OF MOTION 



Consider the Lagrangian density for a massive Abelian vector field plus a scalar field together with Einstein gravity 

C = -^i? + - V{4>) ~ -JF^^F^'^ + -m'A^A^, (3) 

where / is the vector field kinetic function, m is the mass of the vector field, the field strength tensor is 
F^i, = dfj^Aij — diyAfj^ and V{(j>) is the scalar potential. The above can be the Lagrangian density of a massive Abelian 
gauge field, in which case / is the gauge kinetic function. However, we need not restrict ourselves to gauge fields only. 
If no gauge symmetry is considered the argument in support of the above Maxwell type kinetic term is that it is one 
of the few (three) choices [s^l which avoids introducing instabilities, such as ghosts [2^. Also, we note here that a 
massive vector field which is not a gauge field is renormalizable only if it is Abelian |35| . 

It was shown in Ref. [l^ that, as inflation homogenises the vector field diAf^ = 0, the temporal component 
Aq = 0.^ Therefore, without loss of generality, the spacial components of the vector field can be lined-up in the 
z— axis Af^ = (0, 0, 0, ^^(t)). We can then assume a Bianchi-I background with a residual isotropy in the plane 
perpendicular to the vector field expectation value 



ds2 = dt^ - e2a(t) 



2a{t) 



(dx2+ dy2) + e-Mt)dz2l ^ (4) 



where a = e" is the isotropic scale factor and a parameterises the amount of anisotropy. A constant a corresponds to 
a fiat FRW Universe with the Hubble rate given hy H = a. One can shift the values of a and cr by a constant factor 
without changing the physics of the system. 

In Ref. [21| we showed that the scaling / cx a~^^'^ and m (x a leads to scale invariant spectra of super horizon 
perturbations for the vector field components, assuming approximately isotropic quasi-de Sitter inflation. Such per- 
turbations can then contribute to the curvature perturbation of the Universe ^. e.g. through the vector curvaton 
mechanism (l6| . The modulation of the kinetic function / and mass m of the vector field comes from some degree of 
freedom which varies during inflation. The most natural choice for such a degree of freedom is of course the inflaton 
field itself, but other choices are also possible. In this paper, we consider the possibility that the kinetic function and 
mass of the vector field are some function of the inflaton field (j>, i.e. f ~ f {(f>) and m ~ m {(j>). If this is so, then, 
under certain conditions, we demonstrate that the scaling solutions for the kinetic function and mass / cx and 
TO cx a, which give rise to scale invariant perturbation spectra, are in fact dynamical attractor solutions. 



^ If m = then we can set Aq = by a gauge choice. 
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The Einstein equations derived from Eq. ([3]) are 



• 2 -2 

a — a = 



a + 3(i = 



a + 3d(T 



3mp 



3mp 



3mp 



3^(0) - - g^'rn^Al 



-g^-^'fAl+g^-^-m'Al 



(5) 
(6) 
(7) 



where .g^^ = — e 2Q(t)+4(T(t)_ rpj^^ scalar field equation of motion is 

4> + 3a^ + V'{(P)+BA = 0, (8) 

where the prime denotes derivative with respect to </> and Ba is a source term due to the modulation of the kinetic 
function and mass of the vector field from this scalar field. We define this backreaction of the vector field to the scalar 
field dynamics as Ba = Ba-/ + BA-rm where 

Ig^'fAl and BA-,n = -g'^mm' Al, (9) 



Ba- 



A-f 



which we call the /-term and m-term backreaction respectively. 
The equation of motion for the vector field is 

A,+ \a + 4.a+'-\A, + —A, = 0, 

and the vector field energy density is given by pA = Pkin + where 



1 



-y^fA^ 



and 



Va 



-y^m'Al 



The equation describing the acceleration of the Universe is 

a .. . o _ . o 1 

a 



a + a~ 



-02 + V{^) - pki, 



(10) 



(11) 



(12) 



Notice how the vector field potential term Va, has no influence on the acceleration of the Universe, therefore only the 
growth of the vector kinetic term can spoil an inflating Universe. 

It will also be useful to define the following energy density and backreaction ratios 

Ba 



n 



PA 



and 



7^e = 



(13) 



where the energy density of the homogeneous scalar field is given by the usual expression = ^cf' -I- V{(f)). 

The backreaction, Eq. ^ is dynamic with respect to the vector field, as it is a function of A^, but not with respect 
to the scalar field, since it is not a function of (j). Thus, the backreaction may be interpreted as to only modify the 
effective slope of the potential, V^^ = V' + Ba seen by the inflaton. 

If we consider gauge symmetry then / is the gauge kinetic function where / ^ 1 /g'^ and g is the gauge coupling. 
Consequently, because we assume canonical normalisation / — > 1 at the end of inflation, we require that the kinetic 
function is always decreasing in time, f{t) < so that the vector field remains weakly coupled. Note, however, that 
the vector field does not need to be a gauge field necessarily. 

Successful particle production of vector fields during inflation requires that the physical field is effectively massless 
at horizon exit, m/^/f <^ H. However, the vector field can become heavy at some point after horizon exit. In this 
case, particle production can be isotropic |2l[. A heavy vector field oscillates and acts like a pressureless isotropic 
fiuid which does not cause anisotropic stress [13|. Therefore, to consider the effects of a non-zero mass for the vector 
field, we also require the mass to be increasing in time, m{t) > 0, otherwise it will always be negligible. 

Due to the above, we notice that the backreaction always has an opposite sign to the potential slope. This is because 
the scalar field rolls down its potential, so if V'{(j)) > then (p < therefore /'(</>) > and m'{4)) < 0, however if 
F'((^) < then (j) > therefore f'{4>) < and m'{(j)) > 0. This ensures that Ba always has an opposite sign to the 
potential slope, ¥'{(/)). We conclude that the backreaction will always reduce the effective potential slope experienced 
by the inflaton and slow down the scalar field as it rolls down its potential. The effective potential slope as seen by 
the scalar field is then given by 



T/'ff = \/' + 6a = y'(l + 7^B). 



(14) 
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III. DYNAMICAL ANALYSIS 

The Einstein equations Eqs. ©-([T]) together with the scalar and vector field equations Eqs. (|8]) and ([T0| form 
a system of coupled non-linear differential equations. Analytical solutions to this system are extremely difficult if 
not impossible to obtain. However we may understand the qualitative behaviour of the solutions by analysing the 
phase-space of the system. Our aim is to identify attractor solutions which may result from the effect of the vector 
field backreaction on the scalar field dynamics. 

Let us first define a new set of dimensionless, expansion normalised variables 

L = — , X = —= , y = —= , z = , s = —= and fj, = - — ■= . (15) 

en vSmpd vSmpa vSmpa vSmpa c«v/ 

The Friedmann constraint, Eq. ([5]) then becomes 

S^+x^-f y^ + 2^ + s^ = 1. (16) 

We can write the Einstein equations and the scalar and vector field equations as a system of first order coupled 
non-linear differential equations 



dS 
da 

4^ = -3x- Xi/ +Tz'^ -Es'^ + x(3Y,'^ + 3x^ + 2z^ + s^) , (18) 
da ' ' 

dy 

da 

_d£ 

da 

ds 

da 
d/i 

da 



3E + 2z^ - 2s^ -I- S (3E^ -t- 3x^ + 2z^ + s^) , (17) 



Xxy + y (3E2 + 3x^ + 2z^ + s^) , (19) 

-2z - 2zS] - Tzx - /i,s + z (3S]2 + 3x'^ + 2z^ + s^) , (20) 

-s + 2.sE + Esx + fiz + s (3^2 + 3x^ + 2z^ + s^) , (21) 

[(S - r) X + (3S]2 + 3x^ + 2z2 + s^)] fi. (22) 



Where we have defined the dimensionless model parameters, which are the only model-dependent terms in the system 
above, as 

A(a) ^ y|mp (^^) , r(a) ^ ^mp (^^) and S(a) ^ V6mp (^^) . (23) 

Following the arguments at the end of Sec. [TTl we see that A and F have the same sign but opposite to x and S. For 
definiteness we can choose A and F to be positive and x and S to be negative without loss of generality. 
It will also be useful to use the slow-roll parameters defined by 

.(^)^^(^)' = 1a^ and ^(0)^„^^^)=y|„,,pA'-f2e. (24) 

However because of the backreaction affecting the scalar potential slope V'{4>), it will also be useful to define the 
slow-roll parameters in the Hamilton- Jacobi formalism 

i^doono H 1 Ch 

e^ = ——^-— = 3T.^ + 3x^ + 2z^ + s^ and 7?h = = en - ttt — • 25 

a2 2HH 2a en 

In the framework of our expansion normalised variables, we find that 

Id/ 1 dm 

--^=2xF and — = xE. (26) 

/ da m da 

The scaling solutions required to generate a scale-invariant vector field spectrum of perturbations, / cx a"** and m (x a 
are therefore respectively given by the solutions 

2 1 
a; = — — and x = —. (27) 
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The full system above, Eqs. (|17|) - (P^ . is 5-dimensional and the phase-space analysis becomes highly complicated. 
In 2 dimensions, from an analytic point of view one can say a lot more about a system of differential equations. For 
systems of higher dimension one cannot do much more than analyse the stationary points and their stability, see 
Ref. [3^. The stationary points (or critical points) of the system are defined by = -g^ = ... = 0, and their stability 
is established by analysing the behaviour of the solutions close to the critical points. 

A. Stability 

To investigate the stability of the critical point solutions we will consider perturbations about the critical points 

X = Xc + u, y = yc + v , Y, = T,c + w , ••• (28) 
and substituting into the system of equations we obtain the linearized system 









V 




V 


w 


= M 


w 


[ ) 




[ ■■ J 



(29) 



where the matrix M = A4{xc,yc,'^c, ■ ■ ■)■ The general solution for the evolution of linear perturbations can be 
written as 



u{a) 


= Mie™!"- 




hM3e"^" + --- 




v{a) 




- f2e"^" 4 


- fae"^" + • • • 


(30) 


w{a) 


= wie™!"- 


f u;2e"^" 


+ loae"'^" + • • • 





where ut, Vi, wt are constants of integration and the index i gives the dimension of the sub-system under consideration. 

The critical point is considered to be perturbatively stable when the real parts of all the eigenvalues m.i of the matrix 
M are negative. If Im(mi) 7^ then the point is a stable spiral, otherwise it is a node. Generally the eigenvalues will 
be a mixture of positive and negative numbers, in this case the critical point is an unstable saddle point. If however 
the real parts of one or more of the eigenvalues is zero then the point is non-hyperbolic whose stability cannot be 
established using the linearisation procedure above, more sophisticated methods such as the center manifold theorem 
has to be used. In this paper we only consider hyperbolic critical points. 

Now it is important to note that we may write down these solutions to the linearised system Eq. pOp . if and only 
if the eigenvalues m.i, which will generally depend on the model parameters Eqs. (|23p . are constant or slowly varying 
compared to the Hubble scale H. This may be achieved if the model parameters are constant or slowly varying 
themselves, or if the eigenvalues depend weakly on the model parameters. The method to establish the stability of 
the critical points relies on this assumption, and therefore the validity of this assumption has to be verified once we 
have obtained the eigenvalues of the matrix A4. 

In general, we find that the local dynamics of the critical points depend on the dimensionless model parameters, 
A, r and S. When small continuous changes in these model parameters result in dramatic changes in the dynamics, 
the critical point is said to undergo a bifurcation. The values of the parameters which result in a bifurcation at the 
critical point can often be located by examining the linearised system. The bifurcations are located at the parameter 
values for which the real part of an eigenvalue is zero. We will identify the bifurcations in our system as these may 
correspond to interesting changes in the cosmological scenarios described by the system. 

IV. THE MASSLESS VECTOR FIELD CASE 

There has been recent wide interest in the cosmological imp lications of a massless U{1) field whose kinetic function 
is modulated by the inflaton [l^ Authors in Ref. |l9| have shown that for a certain kinetic functional form 
and potential an attractor solution exist which corresponds to a slow-roll phase of inflation that maintains prolonged 
anisotropy. We will show however using dynamical analysis that these scaling solutions can be obtained for a wide 
class of scalar potentials and kinetic functions. 
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Consider a massless vector field, s,/^, 5 = 0. The full system of equations Eqs. PT|) - (P^ . is then reduced to the 
S-dimensional sub-system 

— = -3E + 2 (1 - - - y2) + s r35]2 + 3a-2 + 2 (1 - - - y"^)] , (31) 
da 

4^ = -3a; - Ay2 + T (l - - - i/) + x ISY,'^ + 3x^ + 2 (l - T.^ - x'^ - y^)] , (32) 
da 

4^ = Xxy + y [3E2 + 3x^ + 2 (l - E^ - - y^)! ^ (-33) 
da ' 

together with the Friedmann constraint E^ + a;^ + + = 1. The critical points (xc, ?/c, Zc, Ec) are found by setting 
4^ = 4^ = 4^ = 0. There are 3 critical points for this system, as shown below. 

cia da da ^ ^ ' 

1. The Standard Slow-Roll Solution, SSTZ. 
The first critical point to be considered is 



A fx-' 



SSU: (xe,yc,^c,E,) = I --, yi- ) , 0, I , (34) 

where the eigenvalues of the matrix Ai for this solution are 

mi,2 = -3 + ^ and = -4 + (A + T) . (35) 

Hence the existence and stability of this solution requires 

A < 3 and A (A + T) < 6. (36) 

This critical point becomes the infiationary standard slow-roll solution when A ^ 3 or equivalently e <C 3. The scalar 
field slowly rolls down its potential which dominates the energy density of the Universe j/c ~ 1- We find that the 
vector field energy density vanishes at this critical point Zc = 0, and consequently so does the anisotropy Ec = 0. 
These are completely diluted by the expansion of space and provides an example of the cosmic no-hair theorem (37| . 
The slow-roll parameters Eqs. and are now 

en = e = -A'^ and '/h = ^7 — e, (37) 



and the scalar equation of state is 



2 2 2 



7^ = i + ^^ = ^2:p^ = 9^ ^^^^ 



The method to establish the stability of the system relies on the fact that the eigenvalues are constant or slowly 
varying. This may be achieved if the model parameters are constant or slowly varying, but this may also be achieved 
in the limits A ^ 3 and A (A + F) <C 6 where the eigenvalues depend weakly on the model parameters. This therefore 
ensures the validity of our method in establishing the stability of the critical points. These limits would also guarantee 
the existence and stability of this critical point. 

For slowly varying model parameters we can identify the bifurcations. The bifurcation value A = 3 determines the 
boundary between the existence and non-existence of the critical point. And passing through the bifurcation value 
A(A -I- F) = 6 the stability of the critical point changes. We will see what happens to the other critical points at these 
bifurcations in the following sections. 

2. The Anisotropic Kination Solution, AK,S . 
The second critical point of this system is 

AK.S : {xc, Uc, Zc, Ec) = (x, 0, 0, ^1 - x'^^ ■ (39) 
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This point turns out to be a non-hyperbolic critical point as one of the eigenvalues of the matrix Ai is zero. Therefore 
we cannot establish its stability with the linearised system. In any case we find that this critical point does not 
correspond to a cosmology of interest. This is because the scalar potential vanishes at the point j/c = and therefore 
cannot correspond to an inflating Universe. We also notice that for the anisotropy Ec to be small we require that 
Xc ~ 1, and therefore the Universe is dominated by the scalar kinetic energy. Otherwise we would have a strongly 
anisotropic Universe, contradicting observations. Therefore we will not pursue the analysis of this critical point any 
further. 



3. The Vector Scaling Solution, VSS. 
The vector scaling critical point is given by 

VSS: {xe,yc,z,,^c) = (-6Ai (A + T) , A1A2 + AF + 6 , A1A2 + AF - 6 , 2Ai (A^ + AF - 6)) , (40) 



where Ai ee (SF^ + 4FA + A^ + 12) ^ and A2 = ^3 (3F2 + 2AF - A2 + 12). 

Calculating the eigenvalues of the matrix M for the vector scaling solution becomes highly complicated, however 
we will see that we are interested in the special case where F ^ A and F ^ 1. In these limits the critical point may 
be approximated by 



f V ^ 2 VAT^ 2(Ar-6)\ 
VSS: (xc,yc,2c,Ec) = I --, 1, , — ^ — I 



(41) 



Plots are shown in Appendix (jA ip of the functional dependence of the real parts of the eigenvalues on the model 
parameters. We can clearly see that in the limits F ^ A and F 1 the real parts of the eigenvalues are approximated 

by 

3 

Re[TOi] ~ —3 and Re[m2,3] ~ — -. (42) 

The real parts of the eigenvalues are negative, therefore the critical point is stable i.e. a late-time attractor. In the 
limits considered the eigenvalues depend weakly on the model parameters. This can be seen as plateaus in the plots 
of Appendix (jA 1[) . This therefore ensures the validity of our method in establishing the stability of the critical points. 
The existence and stability of this solution requires 

AF > 6 , F > A and F > 1, (43) 

which are sufficient but not necessary conditions. We can see that the energy density is dominated by the scalar 
potential j/c ~ 1 with the scalar field in slow-roll, therefore we have an inflating Universe. We note that as F — > 00 
the Universe becomes de-Sitter. The energy density ratio in Eq. is given by 

z2 AF-6 

n = ^— ^ ~ « 1. (44) 

The energy density of the vector field tracks that of the scalar field if the dimensionless model parameters are constant. 
Otherwise the energy density ratio is varying, with time-dependence determined by that of the dimensionless model 
parameters. Under the limits considered, the vector field energy density contribution is kept subdominant. The 
backreaction, from Eqs. ([13]) and (|14p leads to the following effective scalar slope 

F /zV 6- AF , 6 



^---aUJ and V:,= j^V'. (45) 



The existence condition AF > 6 from Eq. (|43)) means that the effective slope seen by the scalar field becomes flatter, 
and so the field slows down as it evolves along its potential. This effect is quantified by the slow-roll parameters 

2A , J 2A V6mp /A' F'\ 

eH^^«l and ^ — + —— i- - - \ . (46) 
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Comparing to the SSTZ slow-roll condition Eq. (|37p . we find £}i{VSS) ~ ^ en (557?.). The slow- roll parameter is 
therefore reduced in the vector scaling solution. 
The anisotropy in the expansion, E is given by 

2 /Ar-6\ 2 2 2_ 

The anisotropy is proportional to the energy density ratio and therefore evolves in the same way. This is a generali- 
sation of the result obtained in Ref. [l^ for any potential and kinetic function that satisfy Eq. (|43|) . We can see that 
this vector scaling solution is a new slow-roll inflationary stage which supports small but prolonged anisotropy which 
would otherwise be completely diluted in the standard slow-roll regime. This prolonged anisotropy may be useful 
in generating statistical anisotropy in the curvature perturbation as described in Ref. fl9l . [20j . The non- vanishing 
anisotropy provides a counter-example to the cosmic no-hair conjecture [37| . 

We can now see how the kinetic function behaves at the critical point. As we have obtained the attractor solution 
in Eq. (|?f|) . the kinetic function scales as 

/att(a) <x e-4«. (48) 

The scalar equation of state is 

2x^ 8 

7^ = l + w^^ o^ — ^l. (49) 

+ 2/ 1 

We notice how the equation of state 7^ > and therefore we are not dealing with a phantom held which would 
otherwise violate the strong energy condition. 

We also note that in the VSS parameter space where AF > 6 and F 3> A one of the eigenvalues of the SSTZ, namely 
1713 of Eq. ([55]) is positive. Therefore the SSTZ is an unstable saddle point as long as the condition A < 3 is also 
satisfied, otherwise the SSTZ would not exist. This really tells us that the SSTZ is unstable because the backreaction 
is growing and would eventually effect the dynamics of the system. If however the backreaction is subdominant at 
some early time, the standard slow-roll attractor could be reached prior to the vector scaling solution. This transient 
period of standard slow-roll inflation could last for a large number of e-folds depending on how subdominant the 
backreaction is. The VSS is however the true late-time attractor. 

For the VSS it is more difficult to identify the bifurcations because the eigenvalues are so complicated. However we 
notice that the bifurcation value A(A + F) = 6, observed at the SSTZ, determines the boundary between the existence 
and non-existence of the VSS. We also notice that the stability of the SSTZ and VSS critical points are interchanged 
at the bifurcation. As we are considering the limit F ;5> A for the VSS, then the bifurcation of interest, through which 
the dynamics of the system changes dramatically is AF — 6. 



A. Time Dependence of the Dimensionless Model Parameters 

One important issue that we must deal with is the fact that generally the model parameters, (A, F, S) are time- 
dependent (A(a), F(q!), E!(a)) and therefore the critical points derived above are also time-dependent. The trivial 
case corresponds of course to constant model parameters, where the critical points are the true late time solutions. 
Otherwise they are known as the instantaneous critical points. If the solutions approach the attractors faster than 
the evolution of the critical points themselves, then it is reasonable to assume that the attractor solutions are reached 
and solutions are dragged along with the critical point. Otherwise the late time solutions are determined by the 
asymptotic behaviour of the dimensionless model parameters. In the vicinity of the stationary points the approach 
to the attractors is determined by the eigenvalues of the matrix Ai. For the critical solutions to be reached, their 
time-dependence must be smaller than the eigenvalue of smallest magnitude. We may write this as 



1 dF, 
Fr da 



^6; 



F' 



mp 



< "11,2,3 



(50) 



where Fc = {xc,yc, Zc, Sc'^c)- Therefore we need to establish which variable has the largest time variation and then 
which eigenvalue has the smallest magnitude. 
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1. The Standard Slow-Roll Solution, SSTZ. 

For the SSTZ we see from Eq. (|34|) that Fc = {xcUc), therefore we have to establish which of these two variables 
has the largest time variation. We find that Eq. leads to 



1 dxr 



da 



y|,,p|A'| and 1^ 



(51) 



The existence condition in Eq. ([55)) therefore tells us that the Xc variable has the largest time dependence. 

Next we need to establish which eigenvalue has the smallest magnitude. If the dimensionless model parameters are 
very slowly varying, then from the existence and stability conditions Eq. ([36]), we can see that in the regime F ^ A, 
and considering slow-roll e ^ 1, that |77ii,2| < ItjsI. If we consider the regime F ;:§> A then we find two possibilities, 
first if I < AF < 6 then |mi_2| > |"t-3| otherwise if AF < § then |toi.2| < 1^31. For the standard slow-roll solution, 
Eq. ([50]) with Fc — Xc leads to the conditions 



ssn 



mp\X'\ = \2e~7j\ < 




|-3 + e| 
1-4 + 2e- 



(52) 



The slow-roll conditions e <C 1 and jj^j <C 1 guarantee that the first inequality above is satisfied. Then all we require 
is that ||AF — 4| > 0{ri), to satisfy the second inequality. This can be readily satisfied as long as we are not too 
close to the bifurcation value AF ^ 6. In fact if ||AF — 4| > 0{1), then both inequalities can be met even if \r]\ ~ 1, 
violating this slow- roll condition. In other words, the SSTZ can also apply to fast- roll inflation, where only one of the 
slow-roll conditions e <C 1 is met. 

If the dimensionless model parameters are varying rapidly, then the validity of our method in establishing the 
stability of the critical points requires the strong conditions A <C 3 and AF <C 6. These conditions ensures that the 
eigenvalues depend weakly on the dimensionless model parameters. In this case the eigenvalues of smallest magnitude 
are clearly |mi.2| — 3. Therefore the inequalities above can be readily satisfied even for fast-roll inflation where \ri\ ~ 1. 



2. The Vector Scaling Solution, VSS. 



For the VSS we see from Eq. ([M]) that Fc = {xcUc, Zc,^c), therefore we have to establish which of these four 
variables has the largest time variation. In the limits F ^ A and F ^ 1 we find that Eq. ([5U| leads to 



1 da;^ 



Xc da 



2^6; 



nip 



r 

f2 



1 A^ 
3F2 



1 dye 
yc da 



= \/6i 



mp 



1 dzc 
Zc da 



2V6', 



mp 



f2 



A^ 

f2 



A'F + AF' 



F (AF - 6) 



and 



1 dS, 



•mp 



2 A^ 
3F2 



2 A 
3f 



F2 



(53) 



Sc da 

Therefore the eigenvalue of smallest magnitude, which from Eq. (j42p is |to2,3| — f in the limits considered, has 
to be larger than all of the expressions above, j^-^^ < (or the strong bound <^ 1) for the attractor to be 

reached. Now if we consider a relatively flat potential w|mp|A'| = |r/ — 2e| < 0(1) and because F ^ 1 then clearly 



mp 



< 1. Then if 



AF 76 6 or A'F w -AF', 
we find that all the expressions above can be reduced to a single (strong) condition 



mp 



F' 






f2 




f,2 



< 1. 



(54) 



(55) 



As expected, if the kinetic function is exponential / oc gf^'t'hnp ^ where c is some constant, then the expression above 
is exactly zero. This is because only exponential functions will give constant dimensionless model parameters. Po- 
tentials and kinetic functions obeying F cx 1/A give A'F = — AF'. These have been considered in Ref. [l^ [131 where 



/ cx exp (cmp^ J ■ 



and satisfy the condition in Eq. ([Mj). 
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V. THE MASSIVE VECTOR FIELD 



We now consider a massive vector field whose mass is also modulated by the infiaton, hence s, /i, S 7^ 0. Therefore 
we now have to consider the full, 5-dimensional system of equations Eqs. (|17p - (P^ . 

As we have seen in Sec. II V 31 for a massless vector field, when the /-term backreaction becomes important, the 
dynamics lead to the attractor solution VSS where /att(Qf) oc e~"'". Such scaling for /, results in the generation of a 
superhorizon spectrum of perturbations for the transverse components of the vector field. However, if the vector field 
is massive, we have to consider also the longitudinal component. As shown in Ref. [2l|, to obtain a scale invariant 
spectrum for the latter we need to have m oc a. Thus, we are also interested in finding out whether the attractor 
solution is such that matt(Q?) (x e" simultaneously. But of course there can only be one solution (t>{a) and therefore 
there must be a connection between f{4>) and m{(f>) to obtain the desired scaling of / and m simultaneously. To this 
end we have to impose the following connection 

2 f m ^ ' 

This condition is the price we have to pay to obtain the desired behaviour for both the mass and the kinetic function. 
The question to be investigated is whether, with the contribution of the mass backreaction taken into account, the 
desired scaling remains an attractor solution. 

For the massive vector field case there are now two terms in the backreaction, see Eq. (j9|). Their relative magnitude 
is given by 



Ba- 



Where we considered the connection Eq. ([56|) in the last equality. 

A natural choice of initial conditions for the vector field can be based on energy equipartition grounds. Energy 
equipartition states that at some early time, for example at the onset of inflation, that pkin ~ Va- In our expansion 
normalised set of variables this gives us zf s|, where the subscript 'i' indicates the early epoch at the onset of 

inflation. Considering the connection Eq. ([55]) we find that the initial backreaction ratio is ^^^^ ~ 2. This suggests 

that the m-term backreaction is important from the onset. Therefore we have to be careful as we cannot neglect 
the m-term backreaction even if we are considering a very light vector field. Thus, assuming initial equipartition of 
energy, we cannot simply extrapolate the analysis from the previous massless case where Ba-vi = 0, to a small but 
non-zero massive case, because initially BA~f ^ BA-m- 

Calculating the stationary points and the eigenvalues of the linearised system becomes highly complicated when 
considering a massive vector field. The problem stems from the mass function /x, whose first derivative, taken from 
Eq. ([221), is 

4^ = [(S - r) X -f (3E2 + Zx^ + 2z^ + s^)] fi. (58) 
da 

Under the condition that F has an opposite sign to both S and x we can see that /i is a monotonically increasing 
function. Asymptotically this function will tend to ^ — >■ 00. The phase space analysis is designed to investigate 
the asymptotic (final) behaviour of systems of differential equations. The problem is that we wish to study the 
qualitative behaviour of solutions for a light vector field, which, since the effects of a non-zero mass are increasing in 
time, corresponds to an intermediate phase in the evolution of the system. Therefore, it is not clear whether we can 
understand the system qualitatively in this intermediate regime where the vector field has a small but non-zero mass. 
However, as shown below, it turns out that we can indeed obtain some useful information if we first consider a very 
light field ^ < 1. 



A. The Light Vector Field Case 



In this section we assume that, if the field is light enough, we may neglect any terms proportional to the mass 
function /x in the system of differential equations. We can then establish the existence of any critical points and their 
stability. If any of the stationary points are stable, i.e. attractors, and the mass is small enough as to not effect the 
dynamics, we should expect the solutions to reach these attractors. We later confirm the validity of this assumption 
numerically for different model examples, (e.g. see Fig.[T]). 
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Figure 1: This numerical plot demonstrates the validity of our approximation in neglecting terms proportional to jj. when 
considering a light vector field. The plot shows a numerical solution for the evolution of the kinetic function scaling j-f^ = 2xr, 
given by Eq. (|26p . with respect to the number of elapsed e- folds A'^, for the particular case where the kinetic function, mass 
and scalar potential are exponential. The numerical solution (in blue-thick line) to the full system of equations, Eqs. (|17|) - (|22[ ). 
is indistinguishable from the numerical solution (in red-thin line) where the terms proportional to jj. have been removed until 
the vector field becomes heavy (A^ > 20 in the plot). The plot is also an example of how the vector scaling solution /att oc e"'*" 
is attained in the light and heavy vector field case. Note that the attractors take only a handful of e-folds to be reached, 
corresponding to the oscillating behaviour for TV ~ 5 (light field) and A'^ ~ 20 (heavy field). More examples are shown in 
Sec.lVni 



The full system of differential equations, Eqs. ([T7 1) - (P^ now reduces to the 4-diniensional sub-system given below 

(59) 
(60) 
(61) 
(62) 



dS 






-3S 


da 


da; 






— 3x 


da 


dz 






[-2 


da 




ds 






[-1 


da 



E 3S^ 



2 2 

X — Z 



3x' 

„2\ 



^z' + s-") , 



X [ZY? + 3x2 _^ + 



2S]-ra;+(3S]2 + 3a:2 + 2z2 



3S^ +3a-^ + 2z 



We may obtain a qualitative feel of how the variables z and s evolve by looking at Eqs. (|6ip and (p^ . These 
variables determine the backreaction and vector field energy density. If the energy density is dominated by the scalar 
potential, i.e. y ~ 1, and considering the model parameter F ^ 1, then Eqs. (j6ip and (j62p may be approximated by 



4^ ~ (-2 - Tx) z 
da 



and — — 

da 



{~l + Ex)s = -{-2-rx)s, 



(63) 



where we used the connection Eq. (j56p in the last equality. First we notice that in this regime the z variable evolves 
twice as fast as the s variable. Hence the /-term backreaction evolves two times faster than the m-term backreaction. 
Because of this difference in evolutionary rate of the backreation terms, there is only a small region of parameter 
space where BA-f ^ BA-m a-t the time where the backreaction becomes important to the dynamics of the system. 
Most of parameter space will be taken by either one or the other backreaction terms dominating. 

Let us consider some early time where the vector field energy density and backreaction, and consequently 
anisotropy are completely subdominant. Then a period of slow-roll (or indeed fast-roll) inflation where en ^ 1 
leads to BA-f growing two times faster than BA-m- Remember from Eq. (|25p that the term in the system above 
3^2 4- 2>x^ + 2^2 + = — ^ = Eh and is therefore very small during inflation. The /-term backreaction is therefore 
very likely to completely dominate the total backreaction and will effect the dynamics of the system before the m-term 
does so. We are effectively considering the same system as that of the massless vector field in which only the /-term 
backreaction existed, see Sec. IIVI 

From the approximations in Eq. (|63p we can directly see that at stationary points where z, s 7^ the solutions 
Xc — —^ and/or Xc—h emerge, see Eq. (P7|) . Therefore we observe that both scaling solutions 



/(«) 



cx e 



-4q 



and 



m(a) cx e" 



(64) 
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are inherent in the equations of motions. But of course we require the connection S to ensure that we obtain 

these two solutions simuhaneously. 

We now make these ideas more concrete by looking for stationary solutions to the full system. There are 5 critical 
points for this system. These critical points include the non-hyperbolic anisotropic kination solution, AJCS discussed 
in Sec. IIV 21 We will not consider this point again. The remaining stationary points for the system of Eqs. (|59p - (p^ 
are given below. 



1. The Standard Slow-Roll Solution, SSTZ. 
The standard slow-roll critical point is given by [c.f. Eq. 




(65) 



The eigenvalues of the matrix Ai are 



mi^2 = -S + \x^, mg = -1 + iA(A- S) , and 7714 = -2 + ^A (A + T) . (66) 

The existence and stability of this solution therefore requires 

A < 3 , A (A + r) < 6 and A (A - S) < 3. (67) 

As seen in the massless vector field case, this critical point corresponds to the standard slow-roll inflationary solution 
when A ^ 3. In this region of parameter space the scalar potential dominates the energy density of the Universe 
?/c ~ 1- The vector field energy density and backreaction completely vanishes at the critical point Zc = Sc = 0, and 
consequently so does any anisotropy = 0- AH the results derived from this critical point in the massless case, 
Eqs. ([57)1 and (|38p apply here. The only difference with the massless case is the stability condition. Considering the 
connection S = — ^F, the existence and the stronger stability condition are now 

A < 3 and A (2A + T) < 6. (68) 

These bounds apply for constant or almost constant dimensionless model parameters. Otherwise we require the 
stronger bounds A ^ 3 and AF ^ 6 to make certain that the eigenvalues depend weakly on the dimensionless model 
parameters. This therefore ensures the validity of our method in establishing the stability of the critical points. In 
these limits the existence and stability of this critical point are guaranteed. 

For time-dependent dimensionless model parameters, the motion of the critical point is given by Eq. (jSip . In 
analogy to Eq. (j52p , the constraint on the time-dependence of the dimensionless model parameter becomes 



ssn-. |2e-77| < <^ Imal =|-l + e+fAF| (69) 




If the slow- roll conditions are met e ^ 1 and \ri\ ^ 1, and |AF — 6| > ©(ry), then the inequalities above will be satisfied. 
This can be achieved as long as we are not too close to the bifurcation value AF 9^; 6. The solutions then are able to 
reach and be dragged along with the critical point. 

If the dimensionless model parameters are varying more rapidly, then the validity of our method in establishing the 
stability of the critical points requires the strong conditions A <C 3 and AF <C 6. These conditions ensures that the 
eigenvalues depend weakly on the dimensionless model parameters. In this case the eigenvalue of smallest magnitude 
is clearly 1 771,3 1 — 1- Therefore the inequalities above are readily satisfied if the slow- roll conditions e ^ 1 and |?7| ^ 1 
are met. 



2. The Vector Scaling Solutions, VSS. 

When considering a massive vector field we find two additional vector scaling solutions to the one found in the 
massless field case. Different vector scaling solutions appear because there are now two terms in the backreaction Eq. 
^ which may effect the dynamics of the system. These critical points arise from the effect of /-term backreaction 
only, m-term backreaction only and a combination of the two. We will analyse these points separately as follows. 
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Vector Scaling Solution 1, VSSi. At this critical point tlie vector potential term and tlierefore the m-term 
backreaction vanishes due to Sc = 0. This is the same critical point as the VSS of the massless case [c.f. Eq. (|40|) ]. 

VSSi : (xe,2/c,2c,Se,Sc) =(-6Ai(A + r), A1A2 v/r2 + AF + 6 , A1A2 v/A2 + AF - 6 , , 2Ai (A2 + AF - 6](^9) 

where Ai = (3r2 + 4rA + A2 + 12)"^ and A2 = ^3 (3r2 + 2Ar - A2 + 12). 

The eigenvalues of the matrix M and thus the stability conditions are however different. Considering the connection 
5 = — in Eq. ((56|) . the eigenvalues are given by 

i2 aa /xS 



mi = -A1A2 , m4 = 9Ai (A2 + AE - 4) and 



1 A A 

"^2,3 = --A1A2 



A2± W-8 ( A3E+ ( 2E2 + ^ ) A2+ (e3--eU-^-?E2 



4 / 2 



(71) 



The limits E 3> A and E ^ 1 ensures that we have an inflating Universe ?/c ~ 1 with small anisotropy Sc <C 1, see Eq. 
(|4ip . The critical point is approximated by 



vcrc f V ^ 2 VAE^ 2^AE_6)\ 

VSSi: (xc,j/c,zc,sc,Sc) = I --, 1, — - — ,0, — ^Y2 — )■ y'^) 

The existence of the critical point requires AE > 6. Therefore we are considering the same conditions as those of the 
massless case in Eq. ()43|) . In these limits the real parts of the eigenvalues are approximated by 

3 3AE — 12 
Re[miJ ~ —3 , Re[r?T,2.3j — ^ 2 ' ^ Re[m4j ~ • (73) 

The real parts of eigenvalues toi,2.3 in these limits are the same as for the massless case, see Eq. (|42p . These are 
negative and lead to solutions being attracted to this point. Having now considered a massive vector field we obtain 
an additional eigenvalue 7714 which is very small but positive. We can see this through the existence condition AE > 6 
that < Re [7714] 1. The critical point is therefore technically an unstable saddle point, however because 7774 is so 
small the critical point is said to be quasi-stable. If a solution approaches this critical point it will remain there for a 
long period of time. At the vicinity of the critical point, solutions evolve very slowly. In these limits the eigenvalues 
depend weakly on the model parameters, see Appendix (|A 2[) . This therefore ensures the validity of our method in 
establishing the stability of the critical points. 

All the results found in Sec. IIV 31 for the energy density ratio TZ, backreaction ratio TZb, effective scalar slope V^'j^, 
slow-roll conditions ejj and 7yH, anisotropic stress E and scalar equation of state 7^, in the massless vector field case 
also apply here, see Eqs. (|44|) - (|47|) and (|49|) respectively. However as we have now considered the connection in Eq. 
([55)l . we obtain the attractor solutions simultaneously, see Eq. (|27p 

/att(a) oc e-*°' and m^tt{a) cx e". (74) 

To reach this critical point the /-term backreaction must dominate the total backreaction at the point where it 
becomes dynamically important z'^ ^ s'^. As described above, this is the most likely scenario because the /-term 
backreaction grows twice as fast as the 771-term backreaction if we consider a prior period of standard slow-roll inflation. 

Eor time-dependent model parameters, as seen in the massless case, the motion of the critical point is given by 
Eq. ([53)1 . Therefore, for solutions to be able to reach and be dragged along with the critical point, the eigenvalue of 
smallest magnitude in Eq. ([75]) must be larger than all of the expressions in Eq. ([55)) . We can see from Eq. ([75)) 
that Re[77i4] <C 1. However, when considering the motion of the critical point we do not need to consider this small 
eigenvalue. The reason why is discussed at the end of Sec. lV A"2) Therefore the eigenvalues of smallest magnitude that 
we need to consider are |Re[777,2.3]| — §• Then, if we consider a relatively fiat scalar potential under the constraints in 
Eq. ([54)) . the conditions of Eq. ([53)) . can be reduced to a single condition, see Eq. ([55)) . 

Vector Scaling Solution 2, VSS 2- At this critical point the vector kinetic term and therefore the /-term 
backreaction vanishes due to Zc = 0. The critical point is given by 



VSS 2 : (xe,2/c,^c,Sc,Sc) = (3A3(S-2A), A3A4 V(S2^SAT6), 0, A^A^^i^^EX^, 2A3 (S A - A2 + 3) j 

(75) 
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where A3 = (3 - 5 S A + 2 + 12) ^ and A4 = SW f S2 - -S A + 4 



The case of mterest occurs m the hmits |5| 3> A and |S| ^ 1, this critical point may then be approximated by 
VSS2 : 



{Xc,yc,Zc,Sc,'Sc) — ( — , 1, 0, — 



3S2 



(76) 



The existence of this solution in these limits requires |AS| > 3. We can see that these limits ensure that the scalar 
potential dominates the energy density i/c ~ 1 and therefore leads to an inflating Universe. However these limits 
do not guarantee stability. In Appendix (|A 3| plots are shown of the functional dependence of the real parts of the 
eigenvalues on the model parameters. The real parts of the eigenvalues of the matrix Ai for this critical point in these 
limits are (7714 is obtained by the analytic expression) 



p r 1 T? [ 1 3 3 (AS + 4) 6Ar - 48 
Re[miJ ~ -3, Re[TO2,3j - Re [7714] ~ ^ = ^ , 



(77) 



where we considered the connection Eq. ([55]) in the last equality. With this connection, the critical point in Eq. ((75|) 
is given by 



VSS2 



(78) 



These limits are now F ^ A and F 1. We then find that the real parts of the first three eigenvalues mi_2,3 
are identical to those of VSS 1 above and VSS of the massless case. These eigenvalues drive solutions towards the 
critical point. However, as in VSS 1 above we find an additional eigenvalue 7714 which is very small in the limits 
considered. There is a window 6 < AF < 8 where Re[7Ti4] < and the critical point is stable. Otherwise for AF > 8 
the critical point is quasi-stable. Solutions that approach this critical point evolve very slowly, and therefore remain 
at the critical point for a long period of time. In these limits the eigenvalues depend weakly on the model parameters, 
seen as plateaus in the plots of Appendix (jA 3|) . This therefore ensures the validity of our method in establishing the 
stability of the critical points. 

Having considered the connection Eq. (j56p we find the following result for the energy density ratio 



2(AF-6) 
f2 



« 1, 



(79) 



which is two times larger than the ratio of the massless and VSSi case given in Eq. (|44)) . The energy density of the 
vector field tracks that of the scalar field if the dimensionless model parameters are constant. Otherwise the energy 
density ratio is varying with time-dependence determined by that of the dimensionless model parameters. Under the 
limits considered the vector field energy density contribution is kept subdominant. 

We notice that this critical point leads to the same backreaction ratio TZ^, effective potential slope V^g, slow-roll 



parameters ejj and rjn and equation of state 7^ as the VSS 1 
and (gni) respectively. 
The anisotropy in the expansion, S is given by 

4 ^AF-6 
3 



above and the VSS of the massless case, see Eqs. 



F2 



IScI « 1. 



(80) 



The anisotropy is proportional to the energy density ratio and therefore evolves in the same way. 

As in the VSS 1 case above, having used the connection in Eq. (j56p . we can simultaneously obtain the desired 
attractor behaviour for the kinetic function and mass given in Eq. (j74p . To reach this critical point the r77-terni 



backreaction must dominate the total backreaction at the point where it becomes dynamically important ^ z^. 



For time-dependent model parameters, the motion of the critical point is given by 



1 dXr 



da 



2V6', 



mp 



f2 



2 A^ 
3f2 



1^ 

Vc da 



= 2V6: 



mp 



1 dzc 
Zr da 



= 



1 dsc 
Sc da 



= 2%/6; 



mp 



F' A' 1 AT -f AF' 
f2^ f2 " 2F(AF-6) 



and 



1 dE, 



Sr da 



= 4V6; 



mp 



r 

f2 



5 A' 1 AT AF' 

3f2 ^ 2F(AF- 6) 



(81) 
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For solutions to be able to reach this critical point and be dragged along with it, the eigenvalue of smallest magnitude 
must be larger than all of the expressions above. As in the VSS i case, the eigenvalue of smallest magnitude which 
drives solutions to the vector scaling solutions, from Eq. (|77p is |m2.3| ~ |. Again, we ignore \m4\. The reason why 
is discussed at the end of Sec. IV A 21 Then if we consider a relatively flat scalar potential under the constraints in 
Eq. ([53]), the conditions in Eq. ([5T|) . can be reduced to a single condition, see Eq. ([55]) . 



Vector Scaling Solution 3, VSS3. At this stationary point the kinetic and potential terms of the vector field 
energy density are comparable. Therefore the /-term and m-term backreaction are both important to the dynamics 
of the system. The critical point is given by 



VSSs 



{xc,yc,z„s„J:c) = (-6A5, As V6r2 + 3^ (2r - 3S) - 3S (r - 2S) + 36 , 
A5 _^ (32 + 6r) A - 3S (r + 2S) - 36 , A^y/dX^ + (-6S + 6r) A - 3r (E + 2S) - 36 , -A5 (-2S + A - F) )(82) 



where A5 = - [2 (S - 2A - F)] ^ . 
Considering the connection in Eq. (|56p . and the limits F 3> A and F ^ 1 we can approximate the critical point as 



VSS3 : {xc,yc,Zc,Sc,T,c) 



2 

T 



VAF-8 VAF-4 
V2F ' f 



A 

3f 



(83) 



In Appendix (|A4|) we show plots of the eigenvalues as functions of the model parameters. We can clearly see that in 
these limits the real parts of the eigenvalues of the matrix Ai are approximated by (m^ is obtained by the analytic 
expression) 



Re[mi] 



Re[TO2,: 



and Re [7714] ~ — 12 



(FA - 8) (FA - 4) 
(3FA-20)F2 



(84) 



Three eigenvalues are found to be identical in all the vector scaling solutions for a massless and light vector field, 
namely mi. 2, 3. These eigenvalues are negative and therefore drive solutions towards the VSS. In the same way as the 
other two vector scaling solutions we find an additional eigenvalue 7714 of very small magnitude. For this critical point 
however we find Re[m4] < 0. Therefore the stationary point is stable and the late-time attractor of the system. All 
solutions will eventually converge to this point. In the limits considered F S> A and F ^ 1, the eigenvalues depend 
weakly on the dimensionless model parameters, as seen by plateaus in the eigenvalue plots of Appendix (|A4p . This 
therefore ensures the validity of our method in establishing the stability of the critical points. The existence and 
stability of this solution requires 



AF > 8, 

which are sufficient but not necessary conditions. 
The energy density ratio is now given by 



n 



F» A 



and 



iAF-S 



F> 1, 



F2 



< 1. 



(85) 



(86) 



The energy density of the vector field tracks that of the scalar field if the dimensionless model parameters are constant. 
Otherwise the energy density ratio is varying with time-dependence determined by that of the dimensionless model 
parameters. Under the limits considered the vector field energy density contribution is kept subdominant. We notice 
again that this critical point leads to the same backreaction ratio TZjs, effective scalar slope V^g, slow-roll parameters 
Eh and ijn and scalar equation of state 7^ as the VSS 2.2 above and the VSS of the massless case, see Eqs. pS)) 
and (|49p respectively. The anisotropic stress for the VSS 3 is however given by 



A 

3f 



IScI « 1. 



(87) 



As in the VSS 1^2 cases above, having used the connection in Eq. (|56|) . we can simultaneously obtain the desired 
attractor behaviour for the kinetic function and mass given in Eq. ([74| . 

For time-dependent model parameters, the motion of the VSS 3 critical point is given by 



1 dxr 



Xc da 



= 2%/6; 



mp 



f2 



4 A^ 
3f2 



1 dye 

yc da 



— V6mp 



A^ 

F2 
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1 dzc 
Zr da 



mp 



r 

f2 



4 

3r2 



1 AT + Ar' 



2 r (Ar - 



1 dsc 
Sr da 



mp 



r 

f2 



4 A' 1 AT + Ar' 

3f2 ^ 2r(Ar-4) 



and 



1 dS, 



da 



f2 



4 A^ 
3f2 



A^ 

Ar 



For solutions to be able to reach this critical point and be dragged along with it, the eigenvalue of smallest magnitude 
must be larger than all of the expressions above. As in the VSS 12 cases above, the appropriate eigenvalue of smallest 
magnitude, which drives solutions to the vector scaling solutions, from Eq. is not |m4| but |m2.3| ~ |. We 

discuss why this is so below. As in previous vector scaling solutions we may approximate the conditions above to a 
single condition if we consider a relatively flat scalar potential under the constraints in Eq. ([M)) . see Eq. ([55]) . 

We have found that all three vector scaling solutions in the limits T ^ A and T ^ 1 , and considering the connection 
in Eq. exhibits the same backreaction ratio TZb, effective scalar potential slope V^ff, slow-roll parameters en and 

rjn and equation of state 7^. We also find that all three vector scaling critical points gives us the attractor behaviour 
desired for the kinetic function and mass of the vector field, from Eq. ([27| 



/att(a) oc 



and 



matt (a) oc e° 



(89) 



To which vector scaling critical point will a solution evolve to depends on the relative magnitude of the /-term and 
m-term backreaction at the point where the backreaction becomes dynamically important. If the /-term backreaction 
dominates the total backreaction at the point where it becomes dynamically important, then solutions will approach 
the VSS 2 critical point. In this case we may approximate the system with ^ s^. We are then dealing with 



effectively the same system as that for the massless case, Eqs. (|3ip - (|55|) for which we know the VSS 1 is an attractor. 

For constant or very slowly varying dimensionless model parameters, the small but positive eigenvalue nii of Eq. 
([75)1 corresponds to flows from the VSSj to VSS 3, where s increases and z decreases. Throughout the flow the Xc 
and Uc solutions are approximately conserved, and therefore so does an inflating Universe with the desired scaling 
solutions for the kinetic function and mass. In the limit T ^ 1 this flow is very slow, it takes many e-folds to deviate 
significantly from the critical point. For constant or very slowly varying dimensionless model parameters the late-time 
attractor solution is the VSS 3 (see Fig. [5]). 

In the same way, if the m-term backreaction dominates the total backreaction at the point where it becomes 
dynamically important <C s^, then solutions will approach the VSS2 critical point. However, as for VSS 2, this 
critical point also has an additional positive eigenvalue m4 of Eq. (|77p which drives solutions very slowly towards 
VSS 3. Eigenvalue 7114 is positive if the dimensionless model parameters do not land in the stable window as discussed 
after Eq. ([75]) . otherwise the VSS 2 is itself stable and VSS 3 does not exist. Again, throughout the flow, the Xc and 
Uc solutions are approximately conserved, and therefore so does an inflating universe with the desired scaling solutions 
for the kinetic function and mass. However, the anisotropy varies along the flow, and it can change from positive to 
negative. 

If a solution is driven near VSS 1^2 then the small but positive eigenvalues, m4 for both points, means that the 
solutions are then slowly attracted towards the VSS 3, which is the late-time attractor (see Fig. 

For varying dimensionless model parameters, the motion of the critical point is generally faster than the flow to 
VSS 3. We have argued that flows to the VSS 3 from VSS 1^2 are characterised by the small eigenvalues m^. However 
the approach to any of the VSS critical points is characterised by the other eigenvalues mi^2.3- As long as the 
time dependence of the critical point is smaller than the approach of the solutions towards the attractor then the 
solutions are able to reach and be dragged along with the critical points for a large number of e-folds. As described in 
Sec. lIV Al in the vicinity of the stationary points the approach of solutions to the attractors is given by the eigenvalues 
of the matrix A4. Therefore for the critical solutions to be reached, their time-dependence must be smaller than the 
eigenvalue of smallest magnitude. We will later confirm these arguments numerically for a number of model examples. 

The eigenvalues of smallest magnitude for VSS 1,2,3 from Eqs. ([75)) . ([77|) and ([M|) respectively, are then 
|Re[m2,3]| — I (We have discounted the eigenvalue 1114 as it only corresponds to flows between the vector scaling 
solutions). These eigenvalues have to be larger than all of the expressions describing the motion of the critical points 
Eqs. 15^, dH]) and dHSl). 

To directly reach the VSS 3 critical point with varying dimensionless model parameters the /-term and m-term 
backreaction must be comparable at the point where they become dynamically important. This corresponds to a 
small area in parameter space as discussed previously. 

We notice that all three vector scaling solutions in the limits F ^ A and F ^ 1 and considering the connection Eq. 
()56p have the same Xc and yc solutions, namely Xc — —y and i/c — 1- This is good news for us as no matter which 
critical point a solutions approaches we obtain an inflating universe with the desired scaling solutions for the kinetic 
function and mass, i.e. the ones which result in scale-invariant spectra for the perturbations of all the components of 
the vector field. 



19 




7. 



Figure 2: This numerical plot demonstrates the flow between the three vector scaling solutions found in the light vector field 
case, for the particular example where the kinetic function, mass and scalar potential are exponential. The plot shows how 
different initial conditions can give rise to solutions approaching different vector scaling critical points, depending on the relative 
magnitudes of z and s. The model parameters are chosen so that the VSSs is the late-time attractor. If a solution does not 
directly fall on this critical point, the small but positive eigenvalues m4 for the V55j,a in Eqs. (|73p and (|77[) . and the small 
but negative eigenvalue 7714 for the VSSs in Eq. (|84|l . drive very weak flow of solutions along the arc (shown in red) towards 
the VSS3. However, throughout the flows, the solution x ~ — 2/r is conserved and therefore so is the scaling solution attractor 
/att oc a~*. The figure also depicts the SSTi. critical point which is an unstable saddle point for the parameter space considered. 



B. The Heavy Vector Field Case 

Once we consider a heavy vector field the dimension of the problem increases because terms in the full system 
Eqs. (|17p -(|22 p proportional to [i cannot be ignored. The dynamical analysis becomes very complicated, however we 
can once again reduce the dimension of the system of differential equations by considering a certain approximation. 
Once the vector field becomes heavy it starts to oscillate rapidly and its equation of motion Eq. PH)) is approximated 
by 

A, + —A, ~ 0, (90) 

whose solution gives 

= (91) 

where "()" denotes average over many oscillations. This leads to (pkin) = (Vk) and therefore (z^) = (s^), see Fig. [S] 
As was shown in Ref. when a vector field becomes heavy it oscillates rapidly and behaves like a pressureless 

isotropic fluid. This is confirmed by analysing the stationary solution determined by = 0, from Eq. (fT7|) 

- 31; + 2 (z^) _ 2 (s^) + S {iY? + 3a-2 + 2 {f) + (s^)) = 0. (92) 

We can clearly see that the solution with = then implies (2^) = (s^) which is satisfied by an oscillating field, 
confirmed numerically in Figs. [3] and Using this result, the full system Eqs. p7)) - (P^ reduces to the following 



20 



3-dimensional system 



4^ = -3E + 3E 
da 



2 2\ 
X 1 



aa 



^^ + x^ + \{l-Y?-x^-V^ 



(93) 
(94) 
(95) 



The Friedmann constraint is now + .t^ + + 2 = 1. There are 4 critical points for this system. These critical 
points include the non-hyperbolic anisotropic kination solution, AJCS discussed in Sec. IIV 21 We will therefore not 
consider this point again. The remaining stationary points (xc, Vc, Zc, ^c) given below. 




0.04: Z 



0.02 _ S 



-0.02 



-0.04 



Figure 3: This numerical plot shows the evolution of the dimensionless variables z and s as the vector field becomes heavy. 
The plot demonstrates the validity of our approximation (z'^) = (s^), used when considering a heavy vector field. 



1. The Standard Slow-Roll Solution, SSTZ. 

The standard slow- roll critical point [c.f. Eq. ((65)) ]. is again one of the stationary solutions. It belongs to the 
sub-set where z, s, S = 0. The eigenvalues of the matrix Ai for the heavy field case are now 

TOi,2 = -3+— and TO3 = -3+-A(2A + F-S). (96) 
o o 

The existence and stability of this solution requires 

A < 3 and A (2A F - S) < 9. (97) 

As previously seen, this critical point corresponds to an inflationary solution when A ^ 3. The results derived from 
this critical point in the massless case, Eqs. (|57)) and (|55)) apply here. The only difference with the massless case are 
the stability conditions. Considering the connection in Eq. (j56p . the existence and stability conditions reduce to 

A < 3 and A Qa + F^ < 6. (98) 

For time-dependent dimensionless model parameters, the motion of the critical point is given by Eq. ([5T|l . In 
analogy to Eq. (j52p , the constraint on the time-dependence of the dimensionless model parameters become 



SSn : \2e-,\ < "^f ^ r:T;' , M.l (99) 




If the slow- roll conditions are met e <C 1 and \r]\ <C 1, and |AF — 6| > 0{ii), then the inequalities above are satisfied. 
This can be achieved as long as we are not too close to the bifurcation value AF ,6 6. Solutions will therefore be able 
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to reach and be dragged along with the critical point. If the dimensionless model parameters are varying more rapidly, 
then the vaUdity of our method in establishing the stability of the critical points requires the strong conditions A <C 3 
and Ar ^ 6. These conditions ensures that the eigenvalues depend weakly on the dimensionless model parameters. 
In this case the eigenvalues have approximately the same magnitude |toi,2,3| — 3. Therefore the inequalities above are 
readily satisfied if the slow-roll conditions e ^ 1 and [tjI <C 1 are met, but they can also be satisfied even if \r]\ ~ 1, 
corresponding then to fast-roll inflation. 



2. The Vector Kination Solution, VKS . 
Let us investigate this new critical point given by 

VKS: (xe,yc,Zc,Ec)= Q(r-S) , 0, ^V^g^O^^^, 0^ , (100) 
whose existence is valid for |r — S| < 3. The eigenvalues of the matrix are 

mi,2 = -^ + i(r-S)^ and ma = ^ + ^ (F - S) + i (F - S)^ (101) 

Considering that the parameters A and F have the same sign and opposite to that of S, as discussed at the end of 
Sec. ini it is clear that is positive. Therefore the critical point is an unstable saddle point. 

For constant or very slowly varying dimensionless model parameters the energy density is dominated by either 
the scalar kinetic term Xc or the oscillating vector field depending on the value of |F — S|. However for varying 
dimensionless model parameters we require the stronger condition |F — Sj ^ 3 to ensure that the eigenvalues have a 
small time-dependence. This therefore guarantees the validity of our method in establishing the stability of the critical 
points. In this case the critical point gives us a vector field dominated energy density where the scalar potential and 
anisotropy vanishes. In any case the point cannot correspond to an inflating Universe Uc = 0, and is therefore of no 
interest to us. 



3. The Vector Scaling Solution, VSS. 



The vector scaling critical point, in the heavy vector field case is given by 



VSS : (xc, yc, zc, Sc) = -3A6 , Ag a/s + (F - S) V 2 A (F - S) , Ag ^2 A2 + A (F - S) - 9, , (102) 



where Ag = (2A + F — ^) ^. We notice that when the vector field becomes heavy and starts oscillating the anisotropy 
vanishes Ec — 0, as expected from Ref. [l^ (see Figs. [Hand [S]). The eigenvalues of the matrix M. are 

mi = — 3A6 (A -f F — S) and 

m2,3 = -^^6|3 (F - S) 3A ± [A (54 - IGA^) (F - S) + (F - sf (45 - IGA^) 

-4A(F-S)^-63A^ + 324]H. (103) 



Let us consider the particular case where we set the connection in Eq. ([55]) . and the limits F ^ A and F ^ 1. The 
vector scaling solution, Eq. (|102p is then approximated by 



VSS: (xe,ye,z„E,)= I --,1,^-^-^^ ,0j. (104) 

We can clearly see from Eq. ()103p , that in the limits and connection considered, the real parts of the eigenvalues 
are approximated by 

3 

Re[mi] ~ —3 and Re[m2,3] — ^2 (105) 
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Therefore the critical point is stable in this asymptotic regime. In these limits the eigenvalues depend weakly on the 
dimensionless model parameters. This therefore ensures the validity of our method in establishing the stability of the 
critical points. The sufficient but not necessary conditions for the existence and stability are 

Ar > 6 , r > A and T > 1. (106) 

In these limits the energy density is dominated by the scalar potential j/c ~ 1) therefore we have an inflating 
Universe. We note that as F — > oo the Universe becomes de-Sitter. The energy density ratio Eq. (fT3| is now given by 



(107) 

The energy density of the vector field tracks that of the scalar field if the dimensionless model parameters are constant. 
Otherwise the energy density ratio is varying with time-dependence determined by that of the dimensionless model 
parameters. Under the limits considered the vector field energy density contribution is kept subdominant. We notice 
that this critical point leads to the same backreaction ratio TZs, effective potential slope V^g, slow-roll parameters 
Eh and t^h and equation of state 7^ as all of the vector scaling solutions of the massless and light vector field cases, 
see Eqs. pS)) . and respectively. We have also obtained the solutions of Eq. ([?f|) . which lead to the desired 
attractor behaviour for the kinetic function and mass, given in Eq. ([74]) for the light vector field case. Therefore the 
desired attractive behaviour given by 

/att(Q!) (X e~^" and matt(a) oc e", (108) 
is conserved as the vector field evolves from being light to heavy (see Figs. 2] and [S]). 




Figure 4: This numerical plot shows an example of the evolution of a solution for a massive vector field. As the backreaction 
is initially subdominant the SSTZ critical point is reached, where the anisotropy E vanishes. Solutions may remain there for 
a long period of time until the backreaction becomes important. Solutions will then move to the VSS where E 7^ and the 
attractor solution x^tt = —2/V is obtained. The vector field may then become heavy. As it does so the anisotropy once again 
vanishes but the solution x^At ~ — 2/r is conserved. This is clearly demonstrated by the plots in Fig. [5] 



For time-dependent dimensionless model parameters, the motion of this critical point is given by 

1 dxr 



Xc da 



2V6; 



mp 



f2 



4 A^ 
3f2 



He da 



A^ 

F2 



1 dzc 
Zr da 



= 2V61 



mp 



f2 



4 A' 1 AT + AF' 
3f2 " 2F(AF - 6) 



and 



1 dE, 



Er da 



= 0. 



(109) 
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For solutions to be able to reach this critical point and be dragged along with it, the eigenvalue of smallest magnitude, 
namely |Re[TO2^3] | ~ |, must be larger than all of the expressions above. Now if we consider a relatively flat scalar 
potential under the constraints in Eq. (|54p . the conditions in Eq. (|f 09p . can be reduced to a single condition, see 
Eq. dSSl). 




Figure 5: These plots show a projection of Fig. [4] in tlie x — y plane. They show in increasing detail how the attractor solution 
Xatt ~ — 2/r is conserved as the solution jumps from the light to the heavy vector scaling solution. 



VI. SUMMARY OF RESULTS 



In this paper we have investigated the qualitative behaviour of solutions to a system of coupled non-linear differential 
equations. This system of equations describe an Abelian vector field whose kinetic function and mass are modulated 
by the scalar field which is driving a period of inflation. The only model dependence of the analysis comes from the 
dimensionless model parameters A, F and S. We have investigated a massless, light and heavy vector field separately. 
Generally we have found two sorts of stationary points of interest in the phase-space. These are the so-called standard 
slow-roll SSTZ and vector scaling solutions VSS. We will summarize our results in the following tables. 



A. The Standard Slow-Roll Solutions, SSTZ. 

Let us first consider the standard slow-roll critical point. This point is observed in the massless, light and heavy 
vector field cases, it is given by 



SSTZ : {xc, Uc Zc, Sc, Sc) — 




(110) 



At this critical point the energy density and backreaction of the vector field vanish Zc ~ Sc = 0- We also notice that 
the anisotropic stress also vanishes, Ec = 0, providing an example of the cosmic no- hair theorem [s^l (see Fig.[4|). The 
slow-roll parameters and scalar equation of state parameter are the same for all cases considered. These are given by 

1 2 2 

en = e = -A^ , m = V - and 7^ = -A^ = -en- (111) 

At the SSTZ critical point the energy density of the Universe is dominated by the scalar potential j/c ~ 1 or equivalently 
3mpH^ ~ V{(j)), if A ^ 3. In this case, all we require is that AF < 6 and A|S| < 3 for stability of the critical point. 
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SSTZ Existence Stability 



Massless A < 3 A < 3 , A(A + T) < 6 



Light A < 3 A < 3 , A(A + r) < 6 , 2A(A - H) < 6 



Heavy A < 3 A < 3 , A(2A + T - H) < 9 

Table I: The Standard Slow-Roll Solutions 



Equivalently we may write these conditions as 



e< 3, 



2 ^' /' . 



and 



V 
'V 



< 1. 



(112) 



When considering varying dimensionless model parameters i.e. time-dependent critical points, an additional condition 
arises. This condition ensures that the approach of solutions to the critical point is faster than their evolution. For 
the SSTZ, the slow-roll condition \r]\ <C 1 ensures this strongly, see Eqs. ([5^ . (jHH) and (|M|) for the massless, light and 
heavy vector field cases respectively. These conditions may even be satisfied if the ?/ slow-roll condition is violated 
\r]\ ~ 1. In other words, the SSTZ can also apply to fast-roll inflation (38| . 



B. The Vector Scaling Solutions, VSS. 



The vector scaling solutions VSS are however more interesting. The critical points are too complicated to write 
again here, c.f. Eqs. (|4ip . ([7^ . (|75|) . ([55]) and p04p for the massless, light {VSS i ,2,3) and heavy vector field cases 
respectively. Note that, in the case of a light vector field, VSS 3 corresponds to both the mass and the kinetic function 
backreacting to the infiaton variation, while VSS 2 corresponds to only the mass backreaction being important and 
VSS 1 corresponds to only the kinetic function backreaction being important. As such the latter is equivalent to the 
VSS point of the massless vector field. Throughout the results described below, we have used the connection from 
Eq. ((56|) . This connection may be written as /'// = —Am! /m. 

When considering the limits T ^ A and T 1, the energy density of the Universe is dominated by the scalar 
potential for all the vector scaling solutions, — 1 or equivalently 3mpH^ ~ V {(/)). These limits correspond to 
inflationary solutions with small but non- vanishing anisotropy Sc ^ 1, providing counter examples to the cosmic 
no- hair theorem (s^l- The anisotropy does however vanish if the vector field becomes heavy (see Fig. 0]). It turns 
out that these limits also guarantee the stability of the critical point. Table [TTl summarises our results for the vector 
scaling solutions 



VSS 


Existence 


Stability 


7^ 




Massless 


AF > 6 


F> A, r> 1 


Ar-6 
r2 


i7^ 


Light: VSSi 


AF > 6 


F> A, r» 1 


Ar-6 


i7^ 


Light: VSS2 


AF > 6 


F> A, r> 1 




-|7^ 


Light: VSS:i 


AF > 8 


F> A, r> 1 


3Ar-16 

2r2 


A 

3r 


Heavy 


AF > 6 


F> A, r> 1 


4 /Ar-6\ 
3\ r2 ) 






Table II: The Vector Scaling Solutions 



It is important to note that the stability conditions for the VSS critical points in Table |TT] correspond only to 
three out of the four eigenvalues of the linearised system. The fourth eigenvalue is tiny and its sign determines flows 
between the VSS critical points (see Fig. [2]). Indeed, as explained at the end of Sec. IV A 21 if AF > 8, the VSS 3 is 
completely stable i.e. all eigenvalues are negative in Eq. (|84|) . and therefore VSS 3 is the absolute late-time attractor, 
whereas VSS 1^2 a-re quasi-stable. They have three negative eigenvalues and one very small but positive eigenvalue, 
see Eqs. ([75)) and ([77)) . However if 6 < AF < 8, then VSS3 does not exist, but VSS2 becomes completely stable, i.e. 
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the late-time attractor. see Eq. ([77| . Therefore, as long as AF > 6, there is always a completely stable VSS critical 
point that solutions will eventually converge to. 

The backreaction, from Eqs. (|T3|) and (fT4|) leads to the same effective slope for all cases considered 



6 — Ar , 6 , , , 

7^e=^^^ and v:s=—V'. (113) 

The existence conditions means that the effective slope seen by the scalar field becomes flatter, and so the field slows 
down as it evolves along its potential. The slow-roll parameters and scalar equation of state parameter are also the 
same for all cases considered, and are given by 

2A 2A V6mp /A' T'\ ^ 8 

ch - ^ , m - Y + — ~ [j^Yj ^<t>- f^- (114) 

All the vector scaling solutions for the massless, light and heavy vector field cases, in the limits considered, have 
the same solution namely Xc = — r/2 (see Figs. S] and [S]). As seen from Eq. (P7|) . this solution together with the 
connection considered Eq. (j56p . lead to the desired attractor behaviour for the kinetic function and mass 

/att(Q!) fx e""^" and matt(Q!) oc e". (115) 

This particular scaling for the kinetic function enables the generation of a scale-invariant perturbation spectrum for 
the transverse component of the vector field. The additional scaling for the mass further enables the generation of 
a scale-invariant perturbation spectrum for the longitudinal component of the vector field, as shown in Ref. (2l| . 
Therefore by imposing this connection, we obtain the scaling of the kinetic function and mass which generates scale 
invariant spectra for the perturbations of the vector field. In general, we have shown that the attractor behaviour is 
not destabilised by the mass backreaction and it can be such that leads to scale invariant spectra for all vector field 
components. 

Except for the exponential functional forms of the scalar potential, kinetic function and mass, where the dimen- 
sionless model parameters are constants, we generally find time-dependent critical points. We then require additional 
conditions which ensure that the approach of solutions to the critical point is faster than their evolution, see Eqs. 
(IMl), dll]), dSHl) and pi?g| for the massless VSS and light VSS i , light VSS 2,3 and heavy VSS respectively. Now if 

we consider a relatively flat inflaton potential y^mp|A'| = {t] — 2e| < 0(1) and if AF 56 6 (or AF 9^; 8 for VSS 3) or 
A'F « — AF', we find that all the conditions can be reduced to a single (strong) condition. 



mp 



J,2 1 



< 1. (116) 



As expected, if the kinetic function is exponential / oc e"^"^/™^, where c is some constant, then the expression above 
is exactly zero. This is because only exponential functions give constant dimensionless model parameters. Potentials 
and kinetic functions obeying F cx 1/A result in A'F = — AF'. These have been considered in Refs. [l^ [l^l where 
/ cx exp (cTTip^ J d(/)) and satisfy the conditions in Eq. (j54p . 



VII. SPECIFIC MODEL EXAMPLES 



We now look at specific model examples choosing the functional forms of the scalar potential and the vector kinetic 
function and mass. We obtain constraints on the dimensionless model parameters and their time-dependence for the 
different attractor solutions of interest, namely the standard slow-roll and the vector scaling solutions. We carry out 
the analysis for the different cases separately, i.e. massless, light and heavy vector field cases. 



A. Constant Dimensionless Model Parameters 



In string theory parameters such as the masses or kinetic functions of vector boson fields are determined by the 
compactification scheme, i.e. their values depend on so-called moduli fields. The moduli are scalar fields which 
parameterize the size and shape of extra dimensions. In that sense they are not fundamental scalar fields, but appear 
so from the viewpoint of the 4-dimensional observer. Typically the dependence of masses and couplings on canonically 
normalised (i.e. with canonical kinetic terms) moduli fields is exponential. And exponential potentials for the moduli 
fields are also reasonable. One such modulus can play the role of the inflaton field. 
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Let us assume the following functional forms for the scalar potential and the vector field kinetic function 

V{c^)^Voe^ and /.e^^'*-*-\ (117) 

where c and q are real positive constants. Our dimensionless model parameters in Eq. (j23[) are given by A = \J^q 

and r = 4y^^, and the slow-roll parameters defined in Eq. are e = ^rj = . 

When we come to consider a massive vector field we need to give a functional form for the mass m{(j)). Considering 
the connection Eq. (|56| . the model parameter S now becomes S — ~2y^^. Equivalently we consider the function 

m(</.) =m,e"^^'^"*'\ (118) 

Constant parameters A and F are only achieved by exponential functions of the form shown in Eqs. ()117p and (jllSp . 
The results for the standard slow- roll solution in this model are given in Table IIIII 
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Table III: Constant dimensionless model parameters 



Clearly we require that c < 1 for the stability conditions to be well defined. We also find the following generic 
results for the standard slow-roll solutions, from Eq. (|llip 

1 2 

CH = e, »yH = -^V and 7^ = -e. (119) 

At the SSTZ critical point the vector field energy density and backreaction vanishes Zc = Sc = 0. The energy density 
of the Universe is dominated by the scalar potential ^ 1, ii q <^ VG, corresponding then to a period of inflation.^ 
In this case the slow- roll conditions are satisfied en = ^ 1- Therefore slow- roll infiation is the late-time attractor. 
The slow-roll conditions guarantee the stability of the critical point as long as c ,6 1. In this model however, once the 
solution reaches the attractor, inflation would never end. This is because the slow-roll parameters are constant and 
the slow-roll conditions would never be violated. One would then have to introduce some kind of hybrid mechanism 
or modification of the potential to end the phase of inflation, i.e. to destabilise the critical point. 

The results for the vector scaling solution in this model are given in Table HVl 
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Table IV: Constant dimensionless model parameters 



^ Inflation with an exponential potential is power-law and it takes place for q < •^2 'y^ < 2/3. One can consider that inflation is 
approximately (quasi)de Sitter only if g <C 1 7^ — > 0. Power-law anisotropic inflation has been recently considered in Ref. |33|| . 
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We also find the following generic results for the vector scaling solutions, from Eq. (|114p 

11 2 

eH = -e, "^^^Yc^ ^^^^ ^"^ ^ 3c2^' ^"^^^^ 

When c > 1 the energy density of the vector field and the backreaction grows. Therefore solutions will be driven 
towards the vector scaling attractors. At the critical point the energy density of the vector field tracks that of the 
scalar field at a constant ratio TZ. However inflation is not spoilt as the vector field energy density remains sub- 
dominant. The sufficient but not necessary stability conditions in Table guarantee that the energy density of the 
Universe is dominated by the scalar potential, i/c ~ 1- 

The backreaction, due to the vector field, becomes proportional to the scalar slope V'{(j)) and a negative contribution. 
The backreaction and the effective slope seen by the scalar field, given by Eq. (|113|) is then 

7^B = and V.'fl- =V' + Ba = -V , (121) 

c c 

hence the effective slope becomes flatter and the fleld slows down when the backreaction takes effect. This effect can 
also be seen from the slow-roll conditions en = ^e = ?7H = li which we observe are smaller than that of the 

SSTZ. The VSS therefore corresponds to a new stage of inflation with non-vanishing anisotropy E. The anisotropy 
becomes constant at the critical point but the stability conditions guarantee that it remains small Ec ^ 1. This is 
key to not spoiling observational requirements. The anisotropy however vanishes if the vector field becomes heavy. 

We also notice that for c > 1 the SSTZ is an unstable saddle point. It is unstable because the backreaction grows 
and will eventually effect the dynamics of the system. Therefore the solution can move from the SSTZ to the VSS 
and in the process the scalar field would slow down as it rolls down its potential. The dynamical evolution of the 
model for a variety of initial conditions and choices of the dimensionless model parameters is shown in Fig. [B) 
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Figure 6: These numerical plots show the evolution of the kinetic function scaling 4-^ — 2xr, given by Eq. (|26p . and 
anisotropy E, with respect to the number of elapsing e-folds A'', for a massive vector field, for the constant dimensionless model 
parameters example. The plots show numerical solutions for three different values of the dimensionless model parameters and 
initial conditions. We observe that the 557?. is obtained prior to the VSS for certain initial conditions. At the 557?. the 
kinetic function takes on a variety of scalings given by different values of 2xr. All solutions move then to the VSS where the 
same attractor 2xr — —4 is obtained, giving /att oc a""* and matt oc a, when considering the connection in Eq. (|56|l . The 
plots clearly show that the same attractor 2xT = —4 is also obtained when the vector field becomes heavy. As expected the 
anisotropy E vanishes at the 557? but not at the VSS if the vector field remains light. The anisotropy does however vanish 
once again when the vector field becomes heavy. We also observe that E can be positive or negative depending on which VSS 
the solution approaches nearest, see Table [III The plots of E also show (in dashed lines) the analytic solutions obtained for the 
VSS attractors, see Table [III We can see how well they agree with the full numerical solutions. 



B. The Chaotic Potential 



The archetypal example for an infiationary potential is the chaotic potential. In Ref. [19| the authors considered 
this potential together with the following kinetic function, 

l/(0) = ^m202 /(0) = /,e^(*'"^'\ (122) 
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The dimensionless model parameters Eq. ([23l) are given by A(0) = \ {jlp'^ ^^^"^ = 2cy| ^^^^ where 
is the mass of the scalar field and c is a real positive constant. The slow- roll parameters defined in Eq. ([^1]) are 

When we come to consider a massive vector field we need to give a functional form for the mass m{(j)). Considering 
the connection Eq. (j56p . the model parameter now becomes S = —c^J^ ("^) ' -E^^^i^^^^'^^ly consider the function 



(123) 



£^(^2_02) 

The results for the standard slow- roll solution in this model are given in Table El 
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Table V: Chaotic Potential 

Clearly we require that c < 1 for the stability conditions to be well defined. We also find the following generic 
results for the standard slow-roll solutions, from Eq. (jllip 



and 
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(124) 



The energy density of the Universe is dominated by the scalar potential j/c ~ 1, if f ' corresponding then 

to a period of inflation. In this case, the slow-roll conditions en <C 1 and r/n = are satisfied. Therefore slow-roll 
inflation is the late-time attractor. It is well known that the above potential satisfies the slow-roll conditions in the 
large field regime. The slow-roll conditions guarantee the stability of the critical point as long as c 9^ 1. The period 
of inflation ends when either the existence or stability conditions are violated, i.e. when the slow-roll condition is 
violated. 

As the dimensionless model parameters are varying in this example, the critical point is therefore moving in phase- 
space, and so we have the additional constraints coming from Eqs. ([52]) . ([69]) and ([99|l for the massless, light and 
heavy vector field cases respectively. Considering the slow-roll limit e <C 1, the conditions that must be satisfied are 
respectively 
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(125) 



These constraints are clearly satisfied if the critical point is stable. Therefore the critical point evolves at a slower 
rate than the approach of the solution to the point. The attractor solution can then be obtained and solutions are 
dragged along with the motion of the critical point. 

The results for the vector scaling solution in this model are given in Table IVTl 

We also find the following generic results for the vector scaling solutions, from Eq. (|114p 
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(126) 



The stability conditions guarantee that the scalar potential dominates the energy density j/c ~ 1- Inflation is not 
spoilt as the vector fleld energy density is sub-dominant, 7?. ^ 1 for all cases. However in this model TZ and Ec are 
not exactly constant but are slowly growing as the scalar field rolls down its potential. The backreaction and the 
effective slope seen by the scalar field, given by Eq. pi3p is then 



nE = - 
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and 



Ba = -v. 

c 



(127) 
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Table VL Chaotic Potential 



The backreaction becomes proportional to the scalar slope ¥'{(/)) and a negative contribution. Therefore the effective 
slope as seen by the scalar field becomes flatter and the field slows down. This effect can be seen from the slow-roll 
conditions en = -^e = ^ 1, which are smaller than that of the SSTZ. The VSS corresponds again to a new stage 
of slow-roll inflation with non-vanishing anisotropy. The anisotropy remains small in the large field regime as to not 
spoil observational requirements. The anisotropy will however vanish if the vector field becomes heavy. We find that 
all the results in the massless vector field case are in agreement with the findings of Ref . . 

We also notice that in the case where c > 1 the SSTZ is an unstable saddle point. As the previous model the vector 
field energy density and backreaction are growing and will eventually effect the scalar dynamics. If the solution moves 
from the SSTZ to the VSS the scalar field would slow down as it rolls down its potential. 

As the dimensionless model parameters are varying in this example, the critical point is therefore moving in phase- 
space, and so we have the additional constraints coming from Eqs. ([55]) . ([5T|l . ([55]) and (|109p for the massless, light 
{VSS 1 ,2,3) and heavy vector field cases respectively. All vector scaling solutions have the same eigenvalue of smallest 
magnitude, namely |Re[TO2.3]| ~ |, where we have neglected the small eigenvalues 771,4 because these correspond to 
flows between the scaling solutions as described at the end of Sec. IV A 21 
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Figure 7: These numerical plots show the evolution of the kinetic function scaling j = 2xr, given by Eq. (|26p . and anisotropy 
E, with respect to the number of elapsing e-folds A'', for a massive vector field, for the chaotic potential model example. The 
plots show numerical solutions for three different values of the dimensionless model parameters and initial conditions. We 
observe that the SSTZ is obtained prior to the VSS for certain initial conditions. At the SSTZ the kinetic function takes on a 
variety of scalings given by different values of 2xr. All solutions move then to the VSS where the same attractor 2xT = —4 is 
obtained, giving /att oc a^* and TTiatt oc a, when considering the connection in Eq. (|56[) . The plots clearly show that the same 
attractor 2xT — —4 is also obtained when the vector field becomes heavy. As expected the anisotropy E vanishes at the 557?. 
but not at the VSS if the vector field remains light. We observe that E is increasing for this example because the dimensionless 
model parameters are time-dependent. The anisotropy does however vanish once again when the vector field becomes heavy. 
We also observe that E can be positive or negative depending on which VSS the solution approaches nearest, see Table HIl The 
plots of E also show (in dashed lines) the analytic solutions obtained for the VSS attractors, see Table HIl We can see how well 
they agree with the full numerical solutions. 



For this model we find A'F + XT' = and from the existence and stability condition described above, we find that 
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the dimensionless model parameters satisfy F' ^ A', hence the constraints are approximated by 




(128) 



The first constraint above apphes to the massless, hght (VSSi^s) and heavy vector field cases. The second condition 
appUes to the Ught {VSS2) critical point. These constraints are easily satisfied through the stability condition F A. 
Therefore the critical point evolves at a slower rate than the approach of the solution to the point. The attractor 
solution can then be obtained and solutions are dragged along with the motion of the critical points. The dynamical 
evolution of the model for a variety of initial conditions and choices of the dimensionless model parameters is shown 
in Fig. El 



C. The Supergravity Inspired Potential 



A natural framework where our models can be realised is in the context of supergravity theories. This is because, in 
such theories the bosonic part of the action is determined by three fundamental functions of the fields of the theory: 
the superpotential W, the Kahler potential K and the gauge kinetic function /. The latter is a holoniorphic function 
of the scalar fields of the theory. Hence, if our vector field is a gauge field of a supergravity theory, then it is natural 
to expect that / = /((/>), i.e. the gauge kinetic function is modulated by the inflaton field. 

Now, the F-term scalar potential in supergravity is given by 

Wm + ^Krn 1 ( Wn + ^K^) - 3^1 ^ (129) 



Vf 



where K"^ " is the inverse of the so-called Kahler metric Kr, 



_ d-'K 



subscripts denote derivatives with respect 



to the corresponding fields (i.e. Wn = and Kn = ■§^) and summation over repeated indexes is assumed. The 
kinetic terms of the scalar fields are given by 



(130) 



The scalar fields are canonically normalised if the Kahler potential is minimal, i.e. if it is of the form K = As 
suming a minimal Kahler potential then the Kahler metric becomes Euclidean {Km'n = 
delta) and the above scalar potential can be written as 
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(131) 



Suppose now that the inflaton field corresponds to a supersymmetric flat direction so that W 7^ VF($) (where (f> here 
denotes the corresponding superfield) . In this case the only dependence of Vf on the inflaton is through the Kahler po- 



tential. Considering also that, for a renormalisable superpotential we expect V^"''-^ > '^j , (where V^"*'^ = 
is the F-term scalar potential in global supersymmetry) we can readily observe that the dependence of Vp on the 



dW 



inflaton is only through the exponential prefactor e^^" 
write 



'p . Since for the inflaton field we have K = <i>(f>* 



F(0) ==Foen™p; , (132) 

where we defined the canonically normalised real inflaton field as </> = |$|/-\/2- Below we employ the above scalar 
potential with two different choices for the functional form of the gauge kinetic function. Since we are considering a 
gauge field, we assume that its mass m is obtained through the Higgs mechanism. To consider that the mass is being 
modulated by the inflaton as well, we assume that there is some connection between the inflaton field (p and the Higgs 
field ip responsible for to. An example of such a possibility is to consider a model of smooth hybrid inflation, where 
the inflaton and the Higgs fields are related as i/i oc 0""'^ during inflation 



1. Model 1 



In our first example we choose a power law kinetic function 

/w^/.(|)'°, 



(133) 
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where c is a real positive constant. The dimensionless model parameters in Eq. (|23|l are given by A(0) = yl (^:^^ 

and r(0) = 4cy^ (^^^ . The slow-roU parameters defined in Eq. ^ are e = i {l^Y and = 1 + {^i^Y ■ 

When we come to consider a massive vector field we need to give a functional form for the mass m{(j)). Considering 
the connection Eq. ([55)) . the dimensionless model parameter now becomes S = —2c^J^(^^^^. Equivalently we 
consider the function 



m{<j>) = J j , (134) 

which means that the Higgs field responsible for the mass of the vector boson is related with the inflaton as -0 oc 0~'^. 
The results for the standard slow- roll solution in this model are given in Table IVIII 
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Table VII: SUGRA Potential, Model 1 



Clearly we require that c < 1 for the stability conditions to be well defined, 
results for the standard slow-roll solutions, from Eq. (jllip 



We also find the following generic 
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(135) 



The energy density of the Universe is dominated by the scalar potential yc ~ 1, if (^Tn^^ ^' corresponding then 
to a period of inflation. In the small field regime, the slow-roll condition en 1 is satisfied, but ryn ~ 1 violating the 
second slow-roll condition. This is expected from the form of the scalar potential and corresponds to the well-known 
77-problem of inflation model-building in the context of supergravity. Therefore we are considering a period of fast-roll 
inflation. The slow- roll condition guarantees the stability of the critical point as long as c 76 1 . 

As the dimensionless model parameters are also varying in this example, the critical point is therefore moving in 
phase-space, and so we have the additional constraints coming from Eqs. ([5^ . and for the massless, light 
and heavy vector field cases respectively. Considering the slow- roll limit e <C 1, even with I77I ^ 1, the conditions that 
must be satisfied are 
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c < - , c < and < 



(136) 



Clearly the constraint on the light vector field case cannot be met as we have set c > 0. In the limit c — > the 
motion of the critical point becomes the same as the approach of solutions to it. Therefore, in the light vector field 
case it is unclear whether solutions can reach the SSTZ critical point, numerical simulations are needed to verify this 
(see Fig. [HI). 

The results for the vector scaling solution in this model are given in Table [VlITl 

We also find the following generic results for the vector scaling solutions, from Eq. (|114p 
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(137) 



As in the previous model we can see that inflation is not spoilt as the vector fleld energy density remains sub- 
dominant in the small field limit. The stability conditions guarantee that the scalar potential dominates the energy 
density. The backreaction due to the vector field becomes proportional to the scalar slope V'{(f>) and a negative con- 
tribution, therefore the effective slope as seen by the scalar field becomes flatter and the field slows down, V^'g = ^V. 
This effect can be seen from the slow-roll conditions en = -^e ^ 1 and t^h = ^V- I'^ tie VSS, inflation is now possible 
for r/H ^ 1 if c is large enough. Therefore even if the SSTZ corresponds to fast-roll inflation, it is possible to obtain 
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Table VIIL SUGRA Potential, Model 1 



a period of slow-roll inflation if the solution reaches the VSS attractor. In that way one can overcome the T^-problem 
which plagues supergravity models of inflation. 

The anisotropy remains small as long as the slow- roll condition is met, and vanishes if the vector field becomes 
heavy. In this model and jSd are not exactly constant but are slowly decreasing as the scalar field rolls down 
its potential. As in previous models, when c > I the SSTZ is an unstable saddle point due to the growing vector 
backreaction. 

As the dimensionless model parameters are also varying in this example, the critical point is therefore moving in 
phase-space, and so we have the additional constraints coming from Eqs. ([53|) . ((8T|) . (|88p and (|I09p for the massless, 
light (VSS 1 ^2,3) and heavy vector field cases respectively. All vector scaling solutions have the same eigenvalue 
of smallest magnitude, namely |Re[TO2,3]| — f, where we have neglected the small eigenvalues 7714 because these 
correspond to fiows between the scaling solutions as described at the end of Sec. IV A 21 
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Figure 8: These numerical plots show the evolution of the kinetic function scahng j = 2xr, given by Eq. (|26|) . and anisotropy 
S, with respect to the number of elapsing e-folds A^, for a massive vector field, for the SUGRA inspired model 1 example. 
The plots show numerical solutions for three different values of the dimensionless model parameters and initial conditions. We 
observe that there is no 557?. solution as the scalar potential is too steep and therefore would correspond to a period of fast-roll 
inflation. Despite the wide range of initial conditions, all solutions move to the VSS where the same attractor 2xr = —4 is 
obtained, giving /att oc a~* and matt oc a, when considering the connection in Eq. (|56|) . The plots clearly show that the same 
attractor 22;r — —4 is also obtained when the vector field becomes heavy. As expected the anisotropy S 7^ at the VSS if the 
vector field remains light. We observe that E is decreasing for this example because the dimensionless model parameters are 
time-dependent. The plots of E also show (in dashed lines) the analytic solutions obtained for the VSS attractors, see Table 
im We can see how well they agree with the full numerical solutions. 

For this model we also find that AT -I- XT' — and from the existence and stability condition described above, we 
find that the dimensionless model parameters satisfy F' ^ A', hence the constraints are approximated by 

2 4 

c > - and c > -. (138) 

3 3 

The first constraint above applies to the massless, light {VSS 2 ^3) and heavy vector field cases. The second condition 
applies to the light {VSS 2) critical point. These constraints are readily satisfied through the existence condition. 
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Therefore the critical point evolves at a slower rate than the approach of the solution to the point. The attractor 
solution can then be obtained and the system is dragged along with the motion of the critical points. The dynamical 
evolution of the model for a variety of initial conditions and choices of the dimensionless model parameters is shown 
in Fig. E 



2. Model 2 



We now consider an exponential kinetic function of the form 



(139) 



where c is a real positive constant. The dimensionless model parameters in Eq. (j23p are given by A(0) — , , ^ 

and r(0) ^ 4cy^ {^) ' ^^'^ slow-roU parameters defined in Eq. ^ are e i {li^) ^^^^ = 1 + (^;;^) ■ 

When we come to consider a massive vector field we need to give a functional form for the mass m{(j)). Considering 
the connection Eq. ([56|). the model parameter now becomes S = — 2ci / 1 ( ) . Equivalently we consider the function 



The results for the standard slow-roll solution in this model are given in Table 
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Table IX; SUGRA Potential, Model 2 



Clearly we require that < 1 for the stability conditions to be well defined. We also find the following generic 
results for the standard slow-roll solutions, from Eq. (jllip 
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As in the previous model, the energy density of the Universe is dominated by the scalar potential if yi^^j ^ 6, 

corresponding then to a period of inflation. This potential satisfies the e <C 1 slow-roll condition in the small field 
regime, this also guarantees the stability of the critical point, however we see that rj ^ 1 violating this slow-roll 
condition. Again, this is the well-known yy-problem of infiation in the context of supergravity. Thus, this turns out to 
be a period of fast-roll inflation. 

As the dimensionless model parameters are also varying in this example, the critical point is therefore moving in 
phase-space, and so we have the additional constraints coming from Eqs. ((52|) . ([69|) and ([99|) for the massless, light 
and heavy vector field cases respectively. Considering the slow- roll limit e ^ 1, even with \r]\ ^ 1, the conditions that 
must be satisfied are 



^'<4' 



— I < and c ( — I < ^. 



(142) 



Clearly the constraint on the light vector field case cannot be met as we have set c (^-f-^ > 0. In the limit c (^-f-^ ~^ 
the motion of the critical point becomes the same as the approach of solutions to it. Therefore, in the light vector 
field case it is unclear whether solutions can reach the SSTZ critical point, numerical simulations are needed to verify 
this (see Fig. 11]). 

The results for the vector scaling solution in this model are given in Table 1X1 
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Table X: SUGRA Potential, Model 2 



We also find the following generic results for the vector scaling solutions, from Eq. (|114p 

^« = Kf ) ' ' ^ i (f ) ""'^ ^ ^ (f ) " ^''^^ 

Once again we can see that inflation is not spoilt as the stability conditions guarantee that the scalar potential 
dominates the energy density. However for this model we also notice that 7?. and \R,b\ are deceasing as the scalar 
field evolves down its potential. The backreaction due to the vector field is a negative contribution to the scalar slope 

Kff = F' + S^ = ^^)F'. (144) 

Therefore the effective slope as seen by the scalar field becomes flatter and the field slows down. The slow-roll 
condition ejj ^ 1 is satisfied as long as the existence condition is met. In the VSS it is now possible for t^h ^ 1 if 
^ ^ As seen previously the slow-roll parameters are smaller than that of the SSTZ. In that way. therefore, one 
can overcome the yy-problem. 

The anisotropy remains small under the stability conditions, but is however decreasing as the scalar field evolves 
down its potential. The anisotropy vanishes as the vector field becomes heavy. 



2xr 




A' 



Figure 9: These numerical plots show the evolution of the kinetic function scaling 4 = 2a;r, given by Eq. (|26|) . and anisotropy 
E, with respect to the number of elapsing e-folds A'^, for a massive vector field, for the SUGRA inspired model 2 example. 
The plots show numerical solutions for three different values of the dimensionless model parameters and initial conditions. We 
observe that there is no 557?. solution as the scalar potential is too steep and therefore would correspond to a period of fast-roll 
inflation. Despite the wide range of initial conditions, all solutions move to the VSS where the same attractor 2xr = —4 is 
obtained, giving /att oc a~* and matt oc a, when considering the connection in Eq. (|56[l. The plots clearly show that the same 
attractor 2a;r — —4 is also obtained when the vector field becomes heavy. As expected the anisotropy E 7^ at the VSS if the 
vector field remains light. We observe that E is decreasing for this example because the dimensionless model parameters are 
time-dependent. The plots of E also show (in dashed lines) the analytic solutions obtained for the VSS attractors, see Table 
im We can see how well they agree with the full numerical solutions. 



35 



This critical point is moving in phase space due to the parameter A = X{4>). Therefore we have the additional 
constraints coming from Eqs. ([5^. ((HT|) . ([55)1 and (|109p for the massless, light {VSS i ^2.3) ^^nd heavy vector field cases 
respectively. All vector scaling solutions have the same eigenvalue of smallest magnitude, namely |Re[m2,3]| — |, 
where we have neglected the small eigenvalues because these correspond to flows between the scaling solutions as 
described at the end of Sec. IV A 21 For this model, the parameter F is constant. Then these constraints are reduced 
to 

(j>i 3c (j>i 3c ^ ^ 

The first constraint above applies to the massless and the light {VSS3) vector field cases. The second condition 
applies to the light (ViS5j 2) ^^nd the heavy vector field cases. If these constraints are satisfied, the critical points 
evolve at a slower rate than the approach of the solution to the point. The attractor solution can then be obtained 
and system is dragged along with the motion of the critical points. The dynamical evolution of the model for a variety 
of initial conditions and choices of the dimensionless model parameters is shown in Fig. [3] 



VIII. CONCLUSIONS 



In this paper we studied the dynamics of a model consisting of a massive Abelian vector field, with a Maxwell-type 
kinetic term and with a kinetic function /(0) and mass m{(f>) modulated by a scalar field 0. This scalar field is driving 
a period of infiation. 

Except for very special cases, analytical solutions to the system of field equations together with Einstein's equations 
are either extremely difficult if not impossible to obtain due to their non-linear nature. However we may understand 
the qualitative behaviour of solutions by phase-space analysis. As we are considering systems of 3 or higher dimension 
we cannot do much more than analyse the stationary points and their stability. We have studied the qualitative 
behaviour of solutions for general scalar potentials V{4>) and functional forms of the vector field kinetic function f{<j)) 
and mass m{(j>). We parametrise the model dependence through the so-called dimensionless model parameters defined 

as A = \j\'nip {jv^^ ^ = \f^^P ('^) ^^^'^ = V^rnp (■^^ • We have obtained the critical points and studied their 
stability for three cases; massless, light and heavy vector field. 

Consider firstly the massless vector field case. This case has already been studied by Ref. [l^ [23| for a particular 
model. However our method is more general and produces their results as a special case. We have found 2 critical 
points of interest in this system. The first critical point corresponds to the standard slow-roll infiation attractor. For 
the solution to be infiationary i.e. 3mpH^ ~ ^(0): the slow-roll condition e ^ 1, must be satisfied. At the critical 
point the energy density of the vector field and the backreaction vanishes. The condition for this to occur is AF < 6. 
We also notice how the anisotropic stress, given by E. vanishes at the critical point S = 0, providing an example of 
the cosmic no-hair theorem. Under the conditions described above, the standard slow-roll critical point is stable i.e. 
a late-time attractor solution. 

Only exponential potentials and kinetic functions give a closed system of equations and therefore exact stationary 
points. However, with more general forms for the potential and the kinetic function, critical points become time- 
dependent. The motion of the critical points in phase-space is given by the time-dependence of the dimensionless model 
parameters. If solutions approach the critical points faster than their evolution then we can expect the attractors 
to be reached and solutions to be dragged along with them. The (strong) condition for this to happen leads to the 
second slow-roll condition \r]\ <i; 1, as long as the bifurcation value is not realised AF 76 6. The condition may even 
be satisfied if the 77 slow-roll condition is marginally violated |?/| ~ 1. In other words, the standard slow-roll attractor 
can also apply to fast-roll inflation. 

The second critical point of interest is the so-called vector scaling solution. This critical point arises due to the 
backreaction effects of the vector field on the dynamics of the scalar field. It therefore requires AF > 6. This condition 
guarantees that the backreaction, and therefore also the vector field energy density, grows. However, before the 
vector field energy density grows so much that it spoils inflation, the backreaction on the dynamics of the scalar 
field intervenes and inhibits further growth. The critical point corresponds to an infiationary solution, if F ^ A and 
F ^ 1. These conditions also guarantee that the critical point is stable i.e. a late-time attractor. The vector field 
backreaction has the effect of flattening the effective potential slope experienced by the inflaton V^ff = -^V . The 
evolution of the inflaton consequently slows down with the slow-roll parameter reduced by 2^ en, compared to the 
standard slow-roll case. The effect of the backreaction generates a new stage of slow-roll inflation where the energy 
density of the vector field is kept subdominant TZ ^ 1, where TZ is the vector to scalar energy density ratio. The 
energy density of the vector field tracks that of the scalar field with constant TZ if the dimensionless model parameters 
are constant, otherwise TZ varies slowly depending on the time-dependence of the dimensionless model parameters. 
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The anisotropic stress E. however does not vanish at the critical point. Indeed, there is a residual anisotropy given 
by E ~ I??.. This provides a counter-example to the cosmic no-hair theorem and has attracted some interest. Such a 
prolonged aniso trop y can be used to generate statistical anisotropy in the curvature perturbation and in gravitational 
waves, see Ref. hJ,M. 

In analogy to the standard slow-roll solution, the motion of the vector scaling critical point leads to a second 
'slow-roll' condition given by mp ^ <C 1. This result is derived when considering a flat potential, i.e. [77 — 2e| < 1 

and as long as the bifurcation value is not approximately realised AF 9^; 6. 

The vector scaling attractor gives the solution f^tt oc e~^" in the limits described above. This particular scaling for 
the kinetic function enables the generation of a scale-invariant transverse component of the vector field perturbation 
spectrum as seen in Ref. [Ill, [29| . 

Now, suppose that the vector field has non-zero mass m. To allow particle production to occur the mass of the 
physical vector field M = raj \fj has to be small (M <C H) when the cosmological scales exit the horizon. If the mass 
is originally negligible, in order for it to have any effect we need to consider that its magnitude is increasing during 
inflation. Thus, it is possible that the field becomes heavy (Af > iJ) by the end of inflation. In this case, when 
considering a massive vector field, we may break the problem down into two regimes; light and heavy field. When 
considering a light vector field, as in the massless case, we find two types of critical points which are of interest. The 
standard slow-roll solution is of course still a critical point of the system. The conditions for the vector field energy 
density and backreaction to vanish are now AF < 6 and |AS| < 3. When considering a massive vector field there are 
two sources of backreaction to the scalar field dynamics. And therefore, depending on which dominates first, there 
are different vector scaling solutions. In fact we find three distinct vector scaling solutions. However considering the 
limits F, |S| 3> A and F, |S| ^ 1, all vector scaling solutions are inflationary with attractor solutions either /att oc oT'^ 
or TT^att oc a. 

To simultaneously attain those attractor solutions we impose the connection /'// = ^Am' /m, i.e. S = — ^F. By 
imposing this connection, we obtain the scaling of the kinetic function and mass which generates a scale invariant 
vector field spectrum. If the vector field perturbations contribute to the curvature perturbation ^ then it can generate 
statistical anisotropy in the spectrum and bispectrum of C, throughout the entire range of the cosmological scales. The 
conditions to obtain an inflationary universe with the above scaling for the kinetic function and mass are minimal, 
namely F ^ 1 and F 3> A. Under these limits we have found that there is a small flow of solutions to one particular 
vector scaling solution, but the attractors for the kinetic function and mass are not affected. In the vector scaling 
solutions the vector field contribution to the density generates anisotropic stress and, if it is non-negligible, it may lead 
to anisotropic inflation, which can produce statistically anisotropic perturbations for the inflaton field itself. They 
can also be used as a source of statistical anisotropy in C, as in the massless vector field case. 

Once the vector field becomes heavy it starts to oscillate rapidly [l^. As expected from Ref. [l^, the anisotropy 
vanishes once the field becomes heavy S = 0. As in the previous vector scaling solutions the conditions to obtain an 
inflationary attractor, are F ^ A, F ^ 1 and AF > 6, when considering the connection S = — ^F described above. 

So, in general, we find that if AF > 6 then the backreaction grows and solutions will fiow to a vector scaling solution 
where /att oc a~^. If we then impose the connection /'// = —4m'/m, we also obtain the solution matt oc a. The 
general conditions for these to be attractor solutions are simply F ^ A and F 3> 1. The above conditions correspond 
to the requirement that the backreaction of the vector field onto the roll of the inflaton is not decreasing in time, so 
that it can eventually affect the inflationary dynamics and drive the system to the vector scaling solution(s). The 
bounds on the dimensionless model parameters F and A allow half the parameter space to lead to the attractor solution 
which generates scale invariant spectra for the perturbations of the vector field components. As demonstrated in our 
examples in Sec. IVIIl the parameter space includes natural values for our dimensionless model parameters. 

The next step is to realise this model in the context of realistic theories beyond the standard model. A promising 
possibility is within supergravity theories. Indeed, our vector field could be a gauge boson of some grand unified 
theory (GUT) in the framework of supergravity. In this case, the gauge kinetic function is a holomorphic function of 
the scalar fields of the theory. Hence it is natural to expect that the rolling of the infiaton would modulate the kinetic 
function. Similarly, the mass will be modulated by a Higgs field, which could well be connected to the inflaton through 
the supersymmetric GUT group. In supergravity, Kahler corrections to the scalar potential result in masses of order 
H during inflation (40| . which implies that the inflaton and the Higgs fields are fast- rolling down the potential slopes, 
causing significant variation in the gauge kinetic function, allowing for our vector scaling solution to be attained. As 
we have shown, after the vector scaling attractor is reached, the backreaction of the vector field slows down the roll 
of the infiaton. This gives rise to a slow-roll phase of inflation, overcoming thereby the infamous 77-problem, which 
plagues inflationary models in supergravity. We have clearly demonstrated this in the last two of our examples in 
Sec. IVIIl This is an added bonus to our setup. 
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Appendix A: Eigenvalues 

The following plots show how the real parts of the eigenvalues of the matrix A4 behave as a function of the 
dimensionless model parameters A, F and S, for the different cases considered. In the massive vector field cases, light 
and heavy, we assume the connection of Eq. (j56p . 



1. Massless Vector Field - Vector Scaling Solution, VSS. 

Here we show numerical solutions to the real parts of the eigenvalues of matrix Ai. These correspond to the VSS 
critical point of the massless vector field case, Eq. (|40p . We show numerical solutions because the analytical solutions 
are far too complicated to reproduce here. However we can clearly see from the plots, that in the limits F A and 
F ^ 1, the real parts of the eigenvalues are approximated by [c.f. Eq. ([H])] 

3 

Re[mi] ~ —3 and Re[m2.3] ~ — -. (Al) 

In the limits considered the eigenvalues depend weakly on the dimensionless model parameters. This can be seen as 
plateaus in the plots. This therefore ensures the validity of our method in establishing the stability of the critical 
point. 




2. Light Vector Field - Vector Scaling Solution, VSSi . 

The eigenvalues of the matrix Ai for the VSS i critical point of the light vector field case, Eq. ([70| . are given by 
Eq. (|7ip . It is clear from the analytical solutions that in the limits F ^ A and F 1 the real parts of the eigenvalues 
are approximated by [c.f. Eq. ([75]) ] 

2 3AF 12 

Re[mi] ~ —3 , Re[m2.3] — ~ 2 ' Re[m4] ~ • (A2) 

In the limits considered the eigenvalues mi^2,3 clearly depend weakly on the dimensionless model parameters. Here 
we are showing a numerical solution to the real part of the eigenvalue TO4 to demonstrate that it also depends weakly 
on the dimensionless model parameters. This can be seen as a plateau in the plot, and therefore ensures the validity 
of our method in establishing the stability of the critical point. 
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Figure 11: This plot shows the eigenvalue m4 of matrix A4 for the VSSi critical point in the light vector field case. 

3. Light Vector Field - Vector Scaling Solution, VSS2- 

Here we show numerical solutions to the real parts of the eigenvalues of matrix A4 which correspond to the VSS2 
critical point of the hght vector field case, Eq. ([75]) . We show numerical solutions because the analytical solutions 
are far too complicated to reproduce here. However we can clearly see from the plots, that in the limits |S| ^ A 
and |S| S> 1, the real parts of the eigeirvalues are approximated by [c.f. Eq. ([77| ] (m4 is obtained by the analytic 
expression) 

T, r 1 o o r 1 3 ^ . . 3 (AS + 4) 6Ar - 48 

Re[mi] ~ -3, Re[TO2,3] - Re[m4] ~ - = :^ . (A3) 

In the limits considered the eigenvalues depend weakly on the dimensionless model parameters. This can be seen as 
plateaus in the plots. This therefore ensures the validity of our method in establishing the stability of the critical 
point. 




Figure 12: These plots show the eigenvalues of matrix M for the VSS2 critical point in the light vector field case. 
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4. Light Vector Field - Vector Scaling Solution, VSS3. 

Here we show numerical solutions to the real parts of the eigenvalues of matrix A4 which correspond to the VSS3 
critical point of the light vector field case, Eq. ([5^ . We show numerical solutions because the analytical solutions are 
far too complicated to reproduce here. However we can clearly see from the plots, that in the limits F ^ A and F ^ 1, 
the real parts of the eigenvalues are approximated by [c.f. Eq. ([M|) ] (7714 is obtained by the analytic expression) 

ReK]^-3, Re[m2,3]c^-^ and Re^] -12 i^j^^^l^^^. (A4) 

In the limits considered the eigenvalues depend weakly on the dimensionless model parameters. This can be seen as 
plateaus in the plots. This therefore ensures the validity of our method in establishing the stability of the critical 
point. 




Figure 13: These plots show the eigenvalues of matrix M for the VSS3 critical point in the light vector field case. 
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